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INTRODUCTION 


Most  military,  commercial,  and  recreational  dives  are  “no-stop  dives,”  known  also  as 
“no-decompression  dives.”  Divers  following  the  depth/bottom-time  instructions  for  no¬ 
stop  diving  can  ascend  directly  to  the  surface  with  a  low  probability  of  contracting 
decompression  sickness  (DCS).  In  contrast,  divers  who  remain  at  a  certain  depth 
longer  than  the  maximal  time  prescribed  by  no-stop  instructions  are  obliged  to  ascend 
to  the  surface  by  using  "stops'1  of  prescribed  times  at  prescribed  depths.  In  1993  the 
U.S.  Navy  changed  the  ascent  rate  from  60  to  30  feet  of  seawater/min  for  the  U.S.  Navy 
Standard  Air  Decompression  Table1  (USN57)  and  in  1999  added  no-stop  limits  for  dives 
at  25  and  30  feet  of  seawater,  gauge  (fswg;  1  fsw  =  3.063  kPa;  33.08  fswg  =  2 
atmospheres  absolute). 

The  USN57  decompression  table  and  tables  used  by  other  navies  are  “deterministic”; 
they  carry  an  underlying  assumption  that  dives  are  either  safe  or  unsafe,  an  assumption 
implying  that  there  is  a  threshold  for  developing  DCS.  Deterministic  models  are 
generated  by  an  intuitive  method  in  which  recorded  outcomes  of  test  dives  and 
impressions  from  operational  dives  are  used  informally.  A  more  modern  alternative  is  to 
use  statistical  methods  to  develop  “probabilistic”  models,2’3  which  assume  that  there  is  a 
graded  response  to  decompression  stress,  so  that  even  ostensibly  safe  dives  carry 
small  but  finite  risks. 

Probabilistic  models  are  generated  by  fitting  dive-outcome  data  to  equations  or 
algorithms,  thus  tying  decompression  instructions  to  objective  facts  about  dive 
outcomes.  Probabilistic  models  allow  a  probability  statement  to  be  assigned  to  any  dive 
profile,  so  distinctions  can  be  made  between  particular  DCS  cases  that  occur  in 
relatively  safe  dives  and  those  that  occur  in  hazardous  dives.  When  speaking  of  a 
probabilistic  model,  we  mean  the  probability  equation  or  algorithm  used  to  fit  the 
“calibration  data”  plus  the  values  of  the  parameters  found  by  the  fit.  Our  calibration  data 
—  information  about  whether  divers  were  stricken  with  DCS  after  experimental  dives  — 
are  from  the  U.S.  Navy  Decompression  Database.4,5 

In  1982  Leitch  and  Barnard6  reviewed  no-stop  instructions  used  by  navies  of  the  United 
States,  Great  Britain,  and  France  and  the  test  dives  that  had  been  instrumental  in 
generating  the  three  sets  of  instructions.7-9  Since  that  time  there  have  been  minor 
adjustments  to  the  tables.10-12  It  is  reasonable  to  assume  that  the  no-stop  instructions 
adopted  by  these  various  navies  are  acceptably  safe;  if  this  were  not  so,  the  navies 
would  have  replaced  their  tables.  However,  it  is  not  clear  what  the  probabilities  of  DCS 
are. 

Our  long-term  goal  is  to  develop  a  set  of  nitrogen-based  decompression  tables  that 
correct  safety  and  capability  deficiencies  of  current  U.S.  Navy  tables.  To  that  end,  we 
developed  a  probabilistic  model  for  standard  air  diving13  that  provides  prescriptions  for 
dives  with  decompression  stops  as  well  as  for  no-stop  dives.  Objectives  of  this  report 
are  to  review  existing  no-stop  dive  instructions  in  the  light  of  probabilistic  and 
deterministic  models  that  we  generate  and  to  estimate  probability  of  contracting  DCS 
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after  no-stop  dives  with  various  bottom  times.  The  results  offer  informed  choices  for 
divers  who  may  be  obliged  to  remain  at  depth  longer  than  stipulated  by  current  no-stop 
instructions. 

Bottom  times  lasting  a  day  or  longer  are  considered  “saturation”  dives;  we  assume  that 
body  tissues  come  close  to  equilibrium  with  the  inert  gas  of  the  breathing  mixture  in  24 
hours.  We  deal  with  saturation  dives  in  another  study.14  No-stop  dives  shallower  than 
40  fswg  with  bottom  times  between  4  hr  and  one  day  (240  to  1 ,440  min)  are  used  in 
ship  husbandry;  they  are  in  the  transition  between  standard  dives  and  saturation  dives. 

A  recent  probabilistic  model15  which  we  will  call  “NMRI  ‘93”  led  to  a  recommendation  to 
shorten  the  bottom  times  for  these  no-stop  subsaturation  dives.  Such  shortening  would 
cause  a  major  decrease  in  the  efficiency  of  U.S.  Navy  operations.  Clarification  of  the 
risks  of  subsaturation  diving  could  also  be  important  in  management  of  rescuers  and 
personnel  being  rescued  from  disabled  submarines.  Submariners  may  be  at  depth  long 
enough  to  become  saturated,  whereas  rescuers,  who  may  be  exposed  to  pressure  for 
shorter  times  than  the  submariners,  are  in  a  subsaturation  status. 

The  U.S.  Naw  Decompression  Database 

The  box  at  the  top  of  Figure  1  illustrates  a  heading  that  summarizes  the  actual  dive 
profile  for  a  test  dive  as  it  is  compiled  in  the  U.S.  Navy  Decompression  Database.4,5 
The  top  line  gives  the  dive  depth  in  fswg;  the  bottom  time  in  minutes;  and  time,  in 

minutes,  for  ascent  to  the  surface 
called  "total  decompression  time" 
(TDT).  The  TDT  is  the  sum  of  times 
at  decompression  stops  and  travel 
time.  Next  (between  the  last  comma 
and  the  end  of  the  line)  is  a  coded 
comment  about  the  profile  and  its 
relation  to  other  profiles  in  the  same 
series. 

In  the  second  line  of  the  heading, 

“1 .000”  signifies  the  gas  breathed 
before  the  dive  (1 .000  is  the  code  for 
air),  and  the  next  “1”  indicates  the 
number  of  person-dives  included  in 
this  particular  profile.  The  “1 .0” 
indicates  that  this  diver  suffered  DCS; 
the  value  here  would  be  zero  if  no 
DCS  had  resulted  and  “0.5”  if  the 
DCS  were  considered  a  “marginal” 
case,  one  having  signs  or  symptoms 
that  are  diagnosed  as  DCS  but  not 
severe  enough  to  require  treatment. 
The  “252.0”  in  the  second  line  of  the 


150.0,  73.2,  143.8,  DIV361  1  REPETS=  1 
1.000,  1,1.0,  252.0,312.0, 


0  50  100  150  200  250 


Figure  1 .  Illustration  of  an  entry  in  the  U.S.  Navy 
Decompression  Database.  The  box  shows  a  sample 
heading  for  a  dive  profile  from  file  EDU1180S.  The 
plot  of  depth  below  sea  level  vs.  time  is  from  the  body 
of  the  tabulation,  which  is  not  shown  here. 
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heading  is  T1 ,  the  last  time  since  the  start  of  the  dive  that  the  diver  was  surely  free  of 
signs  or  symptoms  of  DCS,  and  the  “312.0”  is  T2,  the  first  time  that  signs  or  symptoms 
of  DCS  were  noted.  Units  of  T 1  and  T2  are  minutes. 

In  a  database  entry,  the  heading  is  followed  by  a  body  of  information  about  the  dive 
profile.  The  times  and  depths  plotted  in  Figure  1  are  from  the  body.  The  dive  is  to  1 50 
fswg,  with  bottom  time  of  73.2  min  and  TDT  of  143.8  min,  including  six  decompression 
stops.  The  observation  period  extends  beyond  the  graph  to  1 ,662  min.  If  the  diver 
breathes  a  gas  mixture  other  than  air  during  the  dive,  codes  in  the  body  of  the  profile 
indicate  the  mixture  and  the  times  and  depths  at  which  the  breathing  of  the  non-air 
mixture  commenced  and  stopped. 


METHODS 

NO-STOP  SPREADSHEET 

From  the  U.S.  Navy  Decompression  Database4,5  we  prepared  a  “parent’  Microsoft  Excel 
spreadsheet  of  single-level,  nonrepetitive  dives  for  which  the  breathing  gases  are 
various  nitrogen-oxygen  mixtures;  air  is  the  breathing  gas  for  72%  of  the  total  person- 
dives.  In  the  Database  a  separate  report  describes  the  23  source  files,  each  of  which 
contains  a  series  of  entries  that  provide  information  about  one  to  86  persons  who 
followed  a  particular  dive  profile. 

To  prepare  the  spreadsheet,  we  carefully  studied  the  details  of  the  time  and  depth 
profiles  of  each  entry,  using  the  heading  of  the  entry  as  a  guide  (see  Figure  1).  We 
deleted  profiles  having  more  than  one  distinct  bottom  depth  or  an  indistinct  series  of 
bottom  depths.  When  the  recorded  information  in  a  profile  indicated  that  the  heading 
was  inaccurate  or  that  there  was  a  small  deviation  from  a  square-wave  exposure  to 
depth,  we  made  appropriate  corrections  so  that  the  corrected  depth,  bottom  time,  and 
TDT  pattern  corresponded  approximately  to  square-wave  behavior.  Delays  at  the 
beginning  and  end  of  the  dives  necessitated  by  far  the  most  corrections.  Data  files 
DC4D,  DC4W,  and  EDU885A  needed  the  most  corrections.  We  made  corrections  in 
311  profiles  representing  1051  person-dives;  there  were150  person-dive  corrections  for 
depth  (average  absolute  change  =  1 .59  fswg),  732  for  bottom  time  (average  absolute 
change  =  3.74  min),  and  474  for  TDT  depth  (average  absolute  change  =  2.49  min). 

1)  We  took  bottom  time  to  be  the  difference  between  the  time  when  the  divers 
left  a  depth  of  3  fswg  or  shallower  and  the  time  when  the  divers  left  the 
bottom  depth. 

2)  All  profiles  of  EDU885A  needed  correction  for  lags  at  7  fswg  at  the  beginning 
of  the  dives,  with  an  average  correction  of  -2.54  min  +  1 .34  min  standard 
deviation  (SD). 

3)  We  took  the  TDT  to  be  the  difference  between  the  time  when  the  divers  left 
the  bottom  depth  and  the  time  when  they  reached  a  depth  of  3  fswg  or 
shallower. 
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4)  We  adjusted  for  irregularities  in  the  depth  so  that  the  area  under  a  graph  of 
depth  vs.  bottom  time  is  approximately  equal  to  the  area  under  an 
uncorrected  graph: 

■  When  the  summary  heading  did  not  account  for  a  delay  near  the  final 
bottom  depth,  we  made  the  bottom  depth  less. 

■  When  the  summary  heading  did  not  account  for  a  slow  descent  to 
depth,  we  shortened  bottom  time  and/or  decreased  bottom  depth. 

■  When  the  summary  heading  did  not  account  for  small  variations  in 
bottom  depth,  we  took  average  depth. 

To  use  dives  in  which  the  breathing  gas  is  a  nitrogen  mixture  with  an  oxygen  fraction 
different  from  that  of  air,  we  calculate  an  “equivalent  air  depth”  ( EAD ,  with  units  of  fswg): 


EAD  = 


rl  -FQ2' 
{  0.79  , 


•(£>  +  33)- 33 


(D 


In  Equation  1 ,  F02  is  inspired  fraction  of  oxygen,  0.79  is  the  fraction  of  N2  in  air,  D  is 
actual  depth  of  the  dive  in  fswg,  and  33  is  the  atmospheric  pressure  at  the  surface 
expressed  in  fswa  (feet  of  seawater,  absolute).  In  the  remainder  of  this  report,  the  word 
“depth”  refers  either  to  the  actual  depth  (if  divers  breathed  air)  or  to  the  EAD  (if  divers 
breathed  a  nitrogen  mixture  with  an  oxygen  percentage  different  from  that  of  air).  When 
a  fixed  partial  pressure  of  02  was  breathed,  we  calculate  the  F02  for  Equation  1  as  FO2 
=  P02/{  1+  D/33),  where  P02  is  the  fixed  partial  pressure  in  atmospheres  absolute. 

We  assume  that  a  probabilistic  model  intended  to  generate  instructions  for  a  particular 
type  of  dive  will  be  best  if  it  is  calibrated  with  data  from  that  type  of  dive.  Accordingly, 
we  prepared  a  second  spreadsheet  comprised  only  of  no-stop  dives,  the  “No-Stop 
spreadsheet.”  To  do  so,  we  calculated  two  items  for  each  entry  in  the  parent 
spreadsheet:  (1 )  an  "expected"  ascent  time  ( Exp  —  divide  the  bottom  depth  by  the 
ascent  rate  that  the  U.S.  Navy  currently  prescribes,  30  fsw/min),  and  (2)  the  ratio  of 
actual  documented  time  for  ascent  to  the  expected  time  for  ascent  ( TDT/Exp ).  Many 
dives  used  an  ascent  rate  of  60  fsw/min,  so  many  no-stop  entries  have  TDT/Exp  near 
0.5.  Our  definition  of  no-stop  dives  is  based  on  rate  of  ascent,  not  on  absence  of  stops 
in  the  profile;  we  exclude  dives  that  have  long  ascents,  even  though  they  have  no 
definite  decompression  stops.  Also,  a  few  accepted  dives  have  short  stops  or  short 
delays  during  ascent.  Some  TDT/Exp  ratios  for  no-stop  dives  are  above  1 .0  because  of 
idiosyncratic  delays  in  ascent  or  because  some  of  the  chambers  used  for  long-duration 
dives  cannot  be  exhausted  at  a  rate  of  30  fsw/min.  In  the  No-Stop  spreadsheet  we 
exclude  dives  having  the  TDT/Exp  ratio  below  0.46  and  above  1 .5,  with  two  exceptions: 
(1)  we  extend  the  TDT/Exp  ratio  to  1 .6  for  a  few  shallow  dives  (between  37  and  41 
fswg)  with  long  bottom  times  (between  302  and  720  min),  and  (2)  we  are  lenient  with 
saturation  dives,  as  described  in  the  paragraph  that  follows.  Note  that  we  would 
exclude  the  sample  dive  shown  in  Figure  1  because  the  diver  spent  more  than  two 
hours  at  decompression  stops  and  thus  has  a  high  TDT/Exp  ratio. 
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The  No-Stop  spreadsheet  for  calibration  of  our  no-stop  models  includes  246  person- 
dives,  12%  of  the  total,  with  bottom  times  of  1 ,440  min  (1  day)  or  longer.  We  include 
these  no-stop  saturation  dives  in  the  spreadsheet  to  compensate  for  a  lack  of 
subsaturation  dives.  Thus,  to  make  recommendations  about  subsaturation  dives  with 
bottom  times  between  a  half  day  and  a  full  day,  we  are  obliged  to  interpolate  between 
saturation  dives  and  standard  air  dives.  The  No-Stop  spreadsheet  contains  181 
saturation  dives  with  air  breathing  (18  DCS  cases,  average  depth  =  24.4  fswg)  and  65 
person-dives  for  which  the  breathing  gas  is  a  fixed  partial  pressure  of  0.4  atmospheres 
(atm)  of  oxygen  in  nitrogen  (4  DCS  cases,  average  depth  23.3  fswg).  Only  20%  of  the 
saturation  dives  meet  our  criterion  for  TDT/Exp  of  1 .5  or  below.  We  accept  saturation 
dives  that  have  TDT/Exp  up  to  21  (TDT  up  to  14  min,  average  TDT/Exp  =  15.1 ,  average 
depth  =  21 .2  fswg).  We  assume  that  these  longer  ascents  do  not  introduce  bias 
because  tissue  susceptible  to  DCS  for  saturation  dives  releases  excess  gas  slowly. 

Table  1  provides  a  summary  of  the  source  files;  files  containing  saturation  dives  appear 
in  rows  16  through  23,  inclusive.  The  totals  at  the  bottom  show  104  DCS  cases  among 
the  2,037  person-dives  (see  Appendix).  After  rereading  the  original  documentation  for 
file  EDU159A,16  we  revised  the  information  in  it;  the  revised  file  is  EDU159AVL.  Row  22 
in  Table  1  is  for  one  no-stop  saturation  dive  to  38  fswg,  an  air-breathing  dive  that 
resulted  in  DCS  (unpublished  observation  on  file  at  NEDU,  Panama  City,  FL  —  “Outline 
of  Results  Incident  to  Dives  (Compressed  Air  Exposures)  Performed  in  the  Pressure 
Chamber  of  the  Experimental  Diving  Unit,  Dec-May  1939-1940”).  The  right-hand 
column  of  Table  1  gives  information  on  the  TDT/Exp  ratio.  The  incidence  columns  (DCS 
cases/100  person-dives)  are  for  observed  incidence  (%  Obs)  and  predicted  incidence 
(%  Pred).  The  Cases  Pred  and  %  Pred  columns  refer  to  results  from  the  statistical 
fitting  of  data  for  the  probabilistic  P-No-Stop  Model;  they  will  be  examined  in  the 
RESULTS  section. 

Some  source  files  are  specialized;  for  example,  profiles  in  row  17  record  saturation 
dives  having  small  ranges  of  depths  and  times,  whereas  row  3  includes  a  large  range  of 
depths  and  times.  Some  files  are  for  mild  and  some  for  hazardous  exposures:  row  1 
records  only  one  DCS  case  in  289  person-dives,  whereas  row  6  records  26  cases  in 
141  person-dives.  In  the  source  files,  123  person-dives  in  the  No-Stop  dataset  are 
designated  as  “marginal”;  for  our  purposes,  we  assigned  these  to  be  no-DCS  cases. 
Marginal  outcomes  have  the  same  bottom  times  and  depths  as  DCS  cases,  with  only  a 
few  exceptions. 

DISTRIBUTIONS  OF  THE  VARIABLES 

Figure  2  shows  distributions  of  pertinent  variables  in  the  No-Stop  spreadsheet.  In 
Figure  2A,  number  of  person-dives  decreases  as  depth  increases,  and  the  inset  shows 
that  incidence  of  DCS  is  unevenly  distributed  through  the  depth  categories,  with  the 
second  and  fourth  having  incidences  that  are  50%  or  25%  as  high  as  those  in  the  other 
two.  In  Figure  2B,  numbers  of  person-dives  are  unevenly  spread  across  the  range  of 
bottom  times,  with  about  40%  of  the  person-dives  lasting  less  than  one  hour.  There  are 
no  data  for  the  120  to  179  min  category;  the  star  in  the  inset  is  for  zero  DCS  cases 
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TABLE  1.  SOURCE  FILES  FOR  THE  NO-STOP  SPREADSHEET 

Depth,  bottom  time,  and  ratio  columns  show  ranges  of  values  within  the  file 


Person  Cases  Cases 


%  TDT/Exp 


ill 

DC4D 

86 

glggll 

5-75 

289 

1 

1.8 

0.3 

0.6 

0.5- 1.4 

2 

DC4W 

41 

57  -  265 

5-60 

69 

3 

.8 

1.2 

0.5 -1.0 

3 

EDU1351NL 

43 

44-189 

5-205 

143 

2 

3.8 

1.4 

wm 

1.2 

4 

EDU159AVL 

5 

34-36 

720 

11 

5 

1.7 

45.5 

15.5 

1.0 

5 

6 

EDU557 

21 

37-297 

5-270 

104 

0 

1.0 

0.0 

1.0 

0.5  - 1 .5 

EDU849LT2 

74 

100-150 

27-60 

141 

26 

28.8 

18.4 

20.4 

1.2 

7 

EDU849S2 

35 

33-40 

720 

60 

13 

13.8 

21.7 

23.0 

1.2 

8 

EDU885A 

16 

60-190 

12-67 

112 

H 

2.9 

3.6 

2.6 

0.8- 1.2 

HI 

KB 

NMR8697 

229 

30-111 

30  -  240 

477 

11 

10.6 

2.3 

2.2 

0.6 -0.8 

NMR97NOD 

9 

40-41 

199 

103 

3 

1.9 

2.9 

1.8 

1.5 

U 

NMRNSW 

43 

39-62 

82  -  364 

86 

5 

3.5 

5.8 

4.1 

0.5 -1.1 

m 

NSM6HR 

16 

32-40 

359  -  360 

47 

3 

2.0 

I 

I 

0.5  - 1 .6 

KB 

PASA 

3 

102 

35-38 

5 

1 

0.3 

20.0 

6.0 

0.7  - 1 .0 

KB 

RNPL52BL 

7 

50-150 

11-90 

87 

0 

0.9 

0.0 

1.0 

0.5 

KB 

RNPLX50 

16 

31  -41 

240  -  720 

57 

5 

9.0 

8.8 

15.8 

1.0 

ASATARE 

20 

ASATDC 

15 

ASATFR85 

3 

ASATNMR 

34 

ASATTNSM 

13 

EDUAS45 

5 

1 

NMR9209 

15 

Totals 

750 

21  -25  2,862-2,882  65 


26  -  33  1,440  -1454  23 


2,160 


23-30 


20-24  I  4,320  -  6,181  I  50 


26-30 


33 


38 


1 ,440  —  2,160  I  12 


20  -  23  4,284  -  4,400  48 


2,037  104 


5.1 

1.6 -7.9 

16.1 

0.9 -3.4 

8.5 

2.0 -2.6 

4.0 

0.5-21.0 

14.7 

1.4- 1.5 

35.0 

1.8 

60.0 

1.0 

IKE 

9.6-18.0 

104 
+  26.5* 


95%  confidence  interval  of  the  prediction 


among  zero  dives.  The  inset  shows  that  DCS  incidence  is  highest  in  the  right-hand 
category,  which  includes  all  the  saturation  dives  and  subsaturation  dives  at  720  min. 
We  define  subsaturation  dives  as  dives  between  240  and  1 ,440  min;  there  are  no 
person-dives  at  all  between  721  min  and  1,440  min,  half  of  the  possible  range.  The 
only  subsaturation  dives  in  the  dataset  are  261  person-dives  with  7  DCS  cases  at 
bottom  times  of  240  min  to  719  min  and  118  person-dives  with  23  DCS  cases  at  720 
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Depth,  fswg 


Bottom  time,  min 


Ratio,  TDT/Exp 

Figure  2.  Frequency  distributions  of  variables  in  the  No-Stop  spreadsheet.  Height  of  columns  ‘ 
shows  total  numbers  of  person-dives;  height  of  black  sections  at  the  top  shows  numbers  of 
DCS  cases.  Insets  show  incidence  of  DCS  in  the  categories  of  the  main  graphs.  A:  depth 
(including  air-equivalent  depth  for  person-dives  in  which  non-air  N2-02  mixtures  are  breathed). 

B:  bottom  time.  C:  ratio  of  observed  TDT  to  expected  TDT. 

min.  The  saturation  dives  are  weighted  toward  high  DCS  incidence  (8.9%  versus  5.1% 
for  the  entire  dataset),  have  a  wide  range  of  bottom  times  (1 ,440  to  6,181  min),  and 
have  a  small  range  of  depths  or  air-equivalent  depths  (from  20  to  38  fswg).  Figure  2C 
shows  uneven  distribution  of  the  TDT/Exp  ratio;  the  category  with  TDT/Exp  ratio  below 
0.8  includes  the  most  person-dives,  and  the  inset  shows  its  DCS  incidence  is  the 
lowest. 

A  practical  probabilistic  model  gives  instructions  for  achieving  a  low  incidence  of  DCS, 
so  it  is  desirable  to  have  many  low-incidence  dives  in  the  dive-outcome  information 
used  to  calibrate  the  model.  According  to  the  probabilistic  P-No-Stop  Model  we  will 
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describe,  estimated  probability  of  decompression  sickness  ( Pdcs )  is  less  than  3%  for 
62%  of  the  total  number  of  person-dives.  The  distribution  of  person-dives  in  this  low- 
Pdcs  subset  of  the  data  is  approximately  the  same  as  that  in  Figure  2.  This  low- 
incidence  range  includes  17  DCS  cases,  most  of  which  are  for  shallow  dives  with  rapid 
ascent. 

We  use  the  No-Stop  spreadsheet  to  develop  both  probabilistic  and  deterministic 
models. 

Major  models: 

•  P-No-Stop  Model  —  probabilistic,  based  on  the  entire  No-Stop  spreadsheet. 

•  D-No-Stop  Model  —  deterministic,  based  on  the  entire  No-Stop  spreadsheet  plus 
a  few  dives  with  short  times  at  decompression  stops. 

Minor  models: 

•  Sat-Onlv  Model  —  probabilistic,  based  on  the  saturation  dives  in  the  No-Stop 
spreadsheet  and  used  only  for  the  Likelihood  Ratio  statistical  test. 

•  No-Sat  Model  —  probabilistic,  based  on  the  entire  No-Stop  spreadsheet  minus 
saturation  dives  and  used  only  for  the  Likelihood  Ratio  statistical  test. 


PROBABILISTIC  NO-STOP  MODELS 

Equation  2  characterizes  probability  in  the  logistic  regression  paradigm: 

PdCS  —  l+e(-LOGIT )  (2) 


Logistic  regression  is  often  used  with  a  simple  LOGIT  function: 


LOGIT  —  b0  +  ■  Xl  +  b2  •  X-,  +  b2  ■  X 3  +  — 

The  X,  are  variables  in  the  data,  and  fr,  are  parameters  to  be  estimated  by  the  statistical 
method.  The  NONLIN  module  of  SYSTAT,17  a  commercial  statistical  program,  allows  us 
to  use  a  complex  formula  for  the  variable  corresponding  to  Xj  in  Equation  3.  After 
exploring  other  possibilities,  we  decided  upon  a  particular  LOGIT  function: 

LOGIT  =  a  +  b  ■  D  •  [(l  -  e{-c'T) )+  d  •  (l  -  e(~  f  T) )]  (4) 
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Variables  in  Equation  4  are  dive  depth  in  fswg  (D)  and  bottom  time  in  min  (7).  Five 
parameters  (a,  b,  c,  d,  and  f )  are  estimated  by  the  statistical  process. 

In  Equation  4  the  exponential  form  for  the  two  terms  containing  T allows  for  the  effect  of 
bottom  time  to  stabilize  after  long  times,  a  result  consistent  with  the  idea  that  Pdcs  does 
not  depend  on  bottom  time  after  tissues  are  saturated  with  gas.  For  a  given  Pdcs,  the 
product  of  depth  and  bracketed  time  functions  is  equal  to  a  constant  over  the  entire 
range  of  no-stop  diving.  Hempleman18  suggested  a  similar  proposition:  that  the  product 
of  depth  and  the  square  root  of  time  is  the  operative  function. 

To  understand  Equation  4,  consider  the  conventional  ideas  that  the  body  is  composed 
of  distinct  compartments  or  tissues  that  have  various  absorption  rates  for  inert  gases 
and  that  the  underlying  cause  of  DCS  is  supersaturation  of  inert  gas  in  a  “critical” 
compartment  during  ascent  from  depth.  For  deep  dives,  compartments  having  short 
halftimes  are  critical,  so  bottom  times  are  short.  Compartments  with  longer  halftimes 
are  critical  in  shallower  dives,  and  compartments  with  the  longest  halftimes  are  critical  in 
saturation  dives.  Equation  4  accounts  for  uptake  of  inert  gas  in  a  hypothetical  pair  of 
compartments.  According  to  the  parameters  to  be  discussed,  one  compartment  has 
halftime  of  21 .5  minutes;  the  halftime  of  the  second  compartment,  420  minutes,  is 
somewhat  longer  than  the  360-minute  halftime  often  used  in  modeling  of  altitude 
decompression.19  The  parameter  estimation  process  indicates  that,  for  a  given  degree 
of  equilibration,  the  DCS  risk  associated  with  the  slow  compartment  is  2.76  times 
greater  than  the  risk  associated  with  the  fast  compartment. 

The  NONLIN  module  of  SYSTAT  provides  an  option,  the  LOSS,  for  characterizing  an 
error  statement  that  differs  from  the  conventional  least-squared-errors  statement  used  in 
statistics.  Equation  5  is  the  LOSS  statement  we  use;  it  is  in  terms  of  a  single  dive- 
outcome  entry: 

LOSS  =  -DIVES  [ DCS  In  (ESTIMATE)  +  noDCS  In  (1-  ESTIMATE)]  (5) 

In  Equation  5,  DIVES  is  the  number  of  person-dives  in  the  group  test  that  makes  up  the 
dive-outcome  entry,  DCS  indicates  whether  subjects  in  the  group  contracted  DCS,  In 
signifies  natural  logarithms,  ESTIMATE  is  the  estimated  Pdcs,  and  noDCS  indicates  that 
no  subjects  in  the  group  had  DCS.  The  DCS  variable  is  either  a  one  or  a  zero,  because 
divers  who  dived  together  but  did  not  contract  DCS  are  listed  together  in  the  same 
profile  ( DCS  =  0,  DIVERS  =  number  of  divers),  but  divers  with  DCS  are  listed  separately 
( DCS  =  1 ,  Divers  =  1 ).  The  LOSS  function  is  the  sum  of  LOSS  statements  given  by 
Equation  5  over  the  entire  calibration  dataset. 

To  begin  estimating  parameters,  the  SYSTAT  program  uses  Equations  2  and  4  with 
tentative,  user-generated  starting  values  of  the  parameters  to  calculate  an  initial  Pdcs 
for  the  first  dive-outcome  data  point.  Then  Equation  5  calculates  the  LOSS  value  from 
that  Pdcs.  The  procedure  is  repeated  for  all  data  points,  and  the  LOSS  is  summed  over 
all  points.  The  summed  LOSS  is  recalculated  iteratively  to  find  a  minimum  LOSS  by  a 
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quasi-Newton  iteration  process  that  uses  the  first  and  second  derivatives  of  the  LOSS 
function  to  find  new  tentative  parameter  values.  The  iterations  continue  until  differences 
between  successive  tentative  parameter  estimates  reach  a  predetermined  small  value. 
The  minimized  positive  LOSS  is  the  same  as  a  maximized  negative  log  likelihood  (LL) 
for  maximum  likelihood  estimation  of  parameters.  In  addition  to  parameter  estimation, 
the  NONLIN  module  computes  asymptotic  standard  errors  ( ASE)  and  the  asymptotic 
correlation  matrix  of  parameters  by  estimating  the  Hessian  (second  derivative)  matrix. 

Three  kinds  of  95%  confidence  intervals  are  of  interest.  (1)  From  the  limited-sized 
samples  of  person-dives,  we  use  the  binomial  theorem  to  calculate  intervals  for  the 
“true”  incidences  of  DCS  in  the  entire  population.  (2)  The  intervals  for  estimated 
parameters  are  given  by  the  NONLIN  module’s  output.  (3)  From  results  of  the  statistical 
analysis  (parameters,  ASEs,  and  correlation  matrix),  we  use  the  equations  described  by 
Ku20  to  calculate  intervals  for  the  estimated  Pdcs]  with  the  Ku  equations,  the  Pdcs  and 
confidence  intervals  can  be  tabulated  for  the  D-T  pair  for  each  entry  in  the  calibration 
dataset  or  for  entries  in  a  set  of  decompression  instructions. 

Chi-square  tests  assess  whether  differences  between  observed  DCS  cases  in  the 
calibration  data  and  predicted  cases  from  the  model  are  greater  than  expected  from 
chance: 

X2  =  I!  (  Co  -  Cp)2/  Cp  (6) 

where  Co  is  cases  observed  in  a  subset  of  the  data  and  Cp  is  cases  predicted  in  the 
same  subset.  Small  values  of  the  calculated  x2  value  indicate  that  the  model  is  a  good 
fit  to  the  calibration  data. 

DETERMINISTIC  NO-STOP  MODEL 

Figure  3  illustrates  our  method  for  generating  a  deterministic  model  and  also  shows  an 
alternative  to  Figure  2  for  viewing  the  distribution  of  depths  and  times  in  the  calibration 
dataset.  In  the  panels  of  Figure  3,  no-stop  dives  fall  in  a  hyperbolic  shape  on  a  plot  of 
bottom  time  versus  depth;  divers  can  stay  at  a  shallow  depth  for  a  long  time,  at  a 
medium  depth  for  an  intermediate  time,  and  deep  for  a  short  time.  Triangles  indicate 
DCS  cases  that  occurred  after  a  dive  to  the  depth/bottom-time  pair  indicated  by  the 
coordinates.  For  DCS  cases  in  divers  who  breathed  a  gas  mixture  with  oxygen  different 
from  the  fraction  of  O2  in  air  (0.21),  the  triangles  are  accompanied  by  special  symbols 
(see  the  caption  to  Figure  3). 

Two  bull’s-eye  symbols  between  150  and  200  fswg  in  Figure  3,  almost  superimposed  on 
each  other,  correspond  to  3  DCS  cases  associated  with  dives  having  decompression 
stops.  These  bull’s-eye  dives  are  in  the  parent  spreadsheet  but  not  in  the  No-Stop 
calibration  data.  Clinical  descriptions  of  them  in  the  source  files  are  limited.4,5  We 
include  them  for  our  deterministic  model  on  the  assumption  that  if  DCS  occurs  with 
dives  having  short  times  at  decompression  stops,  DCS  is  also  likely  with  no-stop  dives 
having  the  same  depths  and  bottom  times. 
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Figure  3.  Plots  of  the  dive-outcome  data.  Triangles  =  DCS  cases  in  one  subject  or,  rarely, 
two  subjects.  Small  circles  =  no  DCS  in  a  number  of  subjects  (subjects/circle  =  1  to  86, 
average  =  3)-  Triangles  with  large  circles  =  DCS  cases  in  which  the  subject  breathed  a 
mixture  with  F02  below  0.21.  Triangles  with  large  squares  =  DCS  cases  in  which  the  subject 
breathed  a  mixture  with  F02  above  0.21 .  Triangles  with  large  X  =  four  DCS  cases  in  which 
the  subject  breathed  a  fixed  P02.  Bull’s-eye  shapes  =  three  selected  DCS  cases  in  which  the 
profile  included  short  decompression  stops.  A:  all  the  data;  most  small  circles  are  obscured 
by  triangles.  B:  short  bottom-time  range.  C:  shallow-depth  range.  D:  illustration  of  our 
method  of  generating  deterministic  no-stop  instructions  by  fashioning  a  curve  below  or  to  the 
left  of  DCS  cases  seen  in  panel  B. 

dives  having  short  times  at  decompression  stops,  DCS  is  also  likely  with  no-stop  dives 
having  the  same  depths  and  bottom  times. 

Consider  any  particular  depth  in  Figure  3:  theoretically,  a  diver  who  remains  at  that 
depth  for  a  short  time  should  be  free  of  DCS,  so  the  symbol  for  the  dive  on  the  Figure  3 
axes  would  be  a  small  circle.  A  diver  who  remains  at  the  depth  too  long  may  contract 
DCS;  if  so,  the  symbol  would  be  a  triangle.  Thus,  circles  should  predominate  at  the 
lower  edge  of  the  hyperbolic  shape,  and  triangles  should  predominate  at  the  upper 
edge.  Such  a  configuration  is  visible  in  Figure  3B,  where  a  region  containing  only  small 
circles  is  below  another  region  where  there  are  triangles. 
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Since  our  models  are  concerned  only  with  no-stop  dives,  our  deterministic  model  can  be 
very  simple.  We  need  only  to  differentiate  between  profiles  that  provoked  DCS  and 
those  that  did  not.  In  contrast,  deterministic  models  for  dives  that  entail  decompression 
stops  during  ascent  must  track  tissue  gases  in  order  to  introduce  a  decompression  stop 
whenever  a  compartment  is  at  a  dangerous  level  of  supersaturation.  The  principle  of 
our  approach  is  to  separate  DCS-free  dives  from  those  that  cause  DCS  on  a  depth 
versus  bottom-time  plot.  We  fashion  a  curve,  by  eye,  to  the  left  of  all  the  DCS  cases 
that  have  been  observed  in  the  experimental  dive  trials.  The  guiding  rules  for  drawing 
the  curve  are  that  it  should  avoid  the  triangles  and  should  curve  smoothly.  The  method 
is  illustrated  in  Figure  3D,  which  shows  the  DCS  cases  from  Figure  3B.  We  call  the 
deterministic  model  defined  by  the  curve  in  Figure  3D  the  “D-No-Stop  Model”:  it  is 
below  or  to  the  left  of  the  bull’s-eye  symbols  and  all  the  DCS  cases  (triangles)  in  the  No- 
Stop  dataset.  Drawing  this  curve  is  clearly  a  subjective  process  that  is  imprecise 
because  of  the  gaps  between  the  data  points.  We  arbitrarily  limit  the  D-No-Stop  Model 
to  depths  of  40  fswg  and  deeper;  bottom  times  are  between  zero  and  137  min. 

Although  a  deterministic  model  is,  by  definition,  not  concerned  with  Pdcs,  we  can  use  a 
probabilistic  point  of  view  to  speculate  about  our  D-No-Stop  Model.  It  takes  148  dives 
with  no  DCS  cases  to  say  with  95%  confidence  that  incidence  is  below  2%  (one-tailed 
binomial  distribution),  so  some  of  the  symbols  for  DCS-free  dives  (small  circles)  in 
Figure  3  may  actually  be  for  dives  that  have  a  fairly  high  probability  for  causing  DCS. 

The  maximum  number  of  dives  per  circle  on  the  graphs  is  only  86,  so  we  cannot  be 
confident  that  Pdcs  is  less  than  2%  for  any  of  the  small  circles.  A  few  of  the  triangles 
can  be  associated  with  Pdcs  less  than  2%,  but  most  probably  correspond  to  higher 
Pdcs.  In  the  face  of  such  uncertainties,  we  are  assuming  that  any  triangle  indicates  that 
the  probability  of  DCS  is  above  an  acceptable  level  when  we  locate  a  curve  to  the  left  of. 
the  triangles. 

In  all  of  the  data  used  here,  seven  DCS  cases  are  noticeably  to  the  left  of  the  trend  of 
the  triangle  DCS  case  symbols: 

•  The  bull’s-eye  symbols  at  the  lower  right  of  Figure  3D  stand  for  three  DCS 
cases. 

•  Two  triangles  inside  large  circles  are  at  equivalent  depths  of  57  and  85  fswg 
for  divers  breathing  F02  lower  than  0.21 ;  one  of  these  stands  for  two  DCS 
cases  [numbers  (2),  (3),  and  (4)  in  Table  2]. 

■  An  additional  DCS  case  is  off  scale  on  Figure  3D;  it  is  for  a  diver  breathing 
F02  higher  than  0.21  [number  (1)  in  Table  2,  a  triangle  inside  a  large  square 
that  is  visible  in  Figures  3A  and  3C,  equivalent  depth  =  30  fswg,  bottom  time  = 
240  min]. 

Table  2  gives  gives  details  of  the  four  DCS  cases  from  the  seven  in  which  divers 
breathed  non-air  mixtures.  All  cases  in  Table  2  are  from  the  same  source  file,  and  the 
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TABLE  2.  NON-AIR  DCS  CASES  WHICH  ARE  TO  THE  LEFT  OF  THE  TREND  OF 
TRIANGLES  IN  FIGURE  3 


Source  file 

Actual 

depth, 

fswg 

Equiv 

depth, 

fswg 

Bottom 

time, 

min 

TDT, 

min 

Ratio  ' 
TDT/Exo 

T2, 

min 

FOj 

(1) 

NMR8697 

50 

30 

240 

1.2 

0.8 

601 

0.4 

(2) 

NMR8697 

46 

57 

60 

1.1 

0.8 

166 

0.1 

(3) 

NMR8697 

46 

57 

60 

1.1 

0.7 

781 

0.1 

(4) 

NMR8697 

71 

85 

30. 

1.5  — 

0.6 _ 

272 

0.1 

Equiv  depth  =  equivalent  air  depth,  calculated  by  Equation  1. 

TDT  =  total  decompression  time;  equals  travel  time  tor  no-stop  diving. 

TDT/Exp  Ratio  =  observed  TDT  divided  by  TDT  expected  with  an  ascent  rate  of  30 
fsw/min. 

T2  =  first  time  the  subject  has  symptoms  of  DCS,  timed  from  the  start  of  the  dive. 
F02  =  fraction  of  02  in  gas  breathed  at  depth. 


subjects  switched  from  the  fixed  fraction  of  O2  to  air  at  30  fswg  during  ascent.  We 
suspect  that  either  these  particular  dives  are  aberrant  or  that  breathing  a  non-air  mixture 
carries  a  greater  DCS  risk  than  breathing  air  does.  In  either  case,  the  non-air  dives  may 
represent  a  different  population  from  air  dives. 

Six  of  the  seven  DCS  cases  that  are  to  the  left  of  the  trend  determine  the  position  of  the 
D-No-Stop  Model  curve  in  Figure  3D;  they  are  cases  identified  by  bull’s-eyes  and  circled 
triangles  in  the  first  two  bullets  above.  Non-air  dives  represent  half  of  the  dives  that 
determine  the  depth/bottom-time  combinations  of  the  D-No-Stop  Model;  non-air  dives 
also  affect  our  probabilistic  models,  but  their  impact  on  the  statistical  fit  is  diluted 
because  the  bulk  of  the  calibration  dives  are  air  dives. 


RESULTS 


PROBABILISTIC  MODEL 

Table  3  presents  parameters  for  the  probabilistic  “P-No-Stop  Model,”  that  are  estimated 
using  the  No-Stop  spreadsheet.  The  LL  value  is  much  better  than  the  LL  for  the  null 
model,  in  which  Pdcs  =  inc,  the  DCS  incidence  for  the  calibration  dataset.  The 
parameters  are  more  than  six  times  greater  than  their  standard  errors  in  all  cases.  The 
correlation  matrix  indicates  that  some  of  the  parameters  are  correlated  with  each  other, 
especially  parameters  c  and  cf.  An  equation  having  three  exponential  terms  and  no 
weighting  parameter  for  any  of  the  terms  gives  the  same  LL  as  the  P-No-Stop  Model, 
but  the  time  constants  for  two  of  the  exponential  terms  are  the  same,  indicating  to  us 
that  Equation  4,  with  two  exponential  terms  and  weighting  of  the  second  by  parameter 
d,  is  preferable  to  having  three  exponential  terms.  An  equation  with  three  exponential 
terms  in  which  two  were  weighted  did  not  improve  the  LL. 
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TABLE  3.  PARAMETERS  FOR  THE  P-NO-STOP  MODEL 

Equation  4:  WG1T  =  -  +  A ' D '  10  “  ('  - 


Model  =  No-Stop,  Equation  4,  LL  =  -323.41  Null  model,  LL  =  -410.68,  inc=  5.1% 


95%  Confidence 

Parameter 

Estimate  ASE*  Param/ASE 

Interval 

a 

-8.837279  0.382609 

-23.10 

-9.59  -  -8.09 

b 

0.068630  0.008662 

7.92 

0.0516-0.0856 

c 

0.032453  0.005125 

6.33 

0.0224  -  0.0425 

d 

2.760433  0.379541 

7.27 

2.015-3.506 

f 

0.001651  0.000182 

9.06 

0.001293-0.002008 

Correlation 

matrix 

a  b  c 

d 

f 

a 

1 .000000 

b 

-0.571482  1.000000 

c 

0.127428  -0.852794  1.000000 

d 

0.175982  -0.879394  0.944642 

1 .000000 

f 

0.081734  -0.597798  0.58453 

0.495313 

1 .000000 

*  ASE  =  asymptotic  standard  errors 


We  used  Equation  4  to  calibrate  two  additional  models:  the  “Sat-Only  Model,” 
calibrated  with  the  saturation  data  alone,  giving  LL  of  -56.44,  and  a  “No-Sat  Model” 
based  on  1 ,741  person-dives  that  remained  in  the  No-Stop  spreadsheet  after  all  the 
saturation  data  were  removed,  giving  LL  of -265.19.  We  performed  a  Likelihood  Ratio 
(LR)  test21-23  with  these  two  minor  probabilistic  models  plus  the  P-No-Stop  Model  to 
learn  whether  it  is  advisable  to  combine  data  from  saturation  and  nonsaturation  no-stop 
dives.  The  LR  test  determines  whether  two  datasets,  when  combined,  are  described  as 
well  as  they  are  by  two  separate  evaluations.  Our  test  compares  the  sum  of  the  LL 
values  of  the  evaluations  using  saturation  data  alone  and  nonsaturation  data  alone  with 
the  LL  value  of  the  evaluation  using  the  combined  data: 


LR  =  2  [  (LL  Sat-Only  +  LL  No-Sat )  ~  LL  No-Stop  ]  (7) 

The  LR  calculated  by  Equation  7  is  compared  with  the  value  of  the  chi-square 
distribution  with  k  +  j-p  degrees  of  freedom,  where  /rand  j are  the  numbers  of 
parameters  in  the  Sat-Only  and  No-Sat  evaluations,  respectively,  and  p  is  the  number  of 
parameters  in  the  No-Stop  evaluation.  The  evaluation  of  Equation  7  is  LR  =  2  {[(-56.44) 
+  (-265.1 9)]  -  (-323.41 )}  =  3.56.  The  chi-square  value  for  2  +  5  -  5  =  2  degrees  of 
freedom  and  95%  confidence  is  5.99.  The  calculated  value  of  LR  of  3.56  is  less  than 
5.99,  the  chi-square  value,  an  indication  that  the  datasets  can  be  combined. 


14 


1 


Predictions  versus  observations:  P-No-Stop  Model 

The  equation  for  the  least-squares  trend  line  for  a  plot  of  the  Cases  Pred  (Y)  versus  the 
Cases  Obs  (X)  columns  of  Table  1  is  Y  =  1 .0425X  -  0.1791 ,  with  square  of  the 
correlation  coefficient  (R2)  of  0.91 .  The  slope  of  the  line  is  nearly  1 .0  and  the  intercept 
nearly  zero,  indications  that  the  P-No-Stop  Model  is  a  reasonably  good  predictor  of  the 
DCS  cases  in  the  source  files. 


Figure  4  displays  Pdcs  isopleths  calculated  from  the  parameters  of  the  P-No-Stop 
Model  along  with  dive-outcome  data  in  Figure  3.  The  trace  for  2%  Pdcs  in  panel  A  has 
several  triangles  below  or  to  the  left  of  it,  but  the  trace  for  1%  Pdcs  in  panel  B  avoids 
most  triangles;  the  exception  is  for  deep  dives  at  the  right  of  panel  B. 


Depth,  fswg 


Figure  4.  Dive-outcome  data  and  probability  isopleths  calculated  by  the  P-No-Stop 
Model.  As  in  Figure  3,  small  circles  =  no  DCS  in  at  least  one  and  usually  several 
subjects,  triangles  =  DCS  cases,  triangles  with  large  circles  =  subjects  who  breathed  a 
low  F02  mixture,  triangles  with  large  squares  =  subjects  who  breathed  a  high  F02 
mixture,  bull’s-eye  shapes  =  selected  DCS  cases  from  profiles  that  include 
decompression  stops.  A:  heavy  curve  represents  Pdcs  of  2.0%.  B:  curves,  from  left  to 
right:  Pdcs  =  1 , 2,  3,  and  20%. 


Figures  5,  6,  and  8  show  correlations  between  predictions  by  the  P-No-Stop  Model  and 
observations  derived  from  subgroups  of  data.  To  produce  the  figures,  we  sorted  the 
entire  No-Stop  dataset,  or  a  subset  of  it,  by  a  variable  and  then  apportioned  the  person- 
dives  into  subdivisions  (bins)  containing  approximately  equal  numbers  of  person-dives. 
For  example,  there  are  approximately  340  person-dives  in  each  of  the  six  bins  in  Figure 
5. 

In  Figure  5,  predicted  DCS  follows  the  observed  DCS  well  throughout  the  ranges  of  the 
variables.  Chi-square  values  for  the  panels  are  well  below  11.1 ,  the  value  for  95% 
confidence  and  5  degrees  of  freedom:  1.1  (panel  A),  2.9  (panel  B),  and  3.9  (panel  D). 
For  statistical  analysis,  it  is  desirable  to  have  even  distribution  of  incidences  over  the 
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Incidence  pred  &  obs  Incidence  pred  &  obs 


Bottom  time,  min 


Figure  5;  Incidences  of  DCS  as  functions  of  a  variable:  observed  (squares  and 
light  traces)  and  predicted  (heavy  traces).  Dashed  traces  show  95% 
confidence  intervals  of  observed  DCS  cases.  A:  X-axis  is  the  average  depth  in 
the  bin.  B:  X-axis  is  the  average  bottom  time  on  a  scale  that  covers  all  the 
data.  C:  short  bottom-time  region  is  amplified.  D:  X-axis  is  the  ratio  of  actual 
TDT  to  the  TDT  expected  with  an  ascent  of  30  fsw/min;  extreme  right-hand 
points  include  saturation  dives,  which  we  accepted  for  the  calibration  dataset 
even  though  they  have  long  TDTs.  E:  low-ratio  region  is  amplified. 
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ranges  of  the  variables.  In  Figure  5A,  predicted  and  observed  DCS  incidences  as 
functions  of  depth  fluctuate  markedly  on  either  side  of  the  overall  5.1%  incidence  in  the 
dataset.  In  Figure  5B,  the  incidences  range  from  less  than  1%  to  more  than  12%.  The 
bottom-time  scale  in  Figure  5B  is  large  enough  to  show  a  bin  that  contains  saturation 
data  at  the  right;  as  a  consequence,  the  other  points  are  squeezed  against  the  Y-axis. 

In  Figure  5C,  the  short-dive  region  is  amplified;  incidence  is  low  for  dives  with  short 
bottom  times.  In  later  sections  we  discuss  evidence  that  deep  no-stop  dives  are  not  as 
safe  as  this  low  incidence  suggests. 

Figures  5D  and  5E  shed  additional  light  on  the  poor  distribution  of  the  TDT/Exp  ratio  in 
Figure  2C.  Incidence  is  very  low  when  the  ratio  is  near  0.5  (for  an  ascent  rate  of  60 
fsw/min)  and  high  when  the  ratio  is  near  1 .0  (for  an  ascent  rate  of  30  fswg/min).  When 
we  divide  the  nonsaturation  person-dives  according  to  TDT/Exp  ratio,  we  find  1,055 
dives  with  the  ratio  between  0.5  and  0.74  (average  TDT/Exp  =  0.6  with  observed 
incidence  of  2.3%)  and  800  dives  with  the  ratio  between  0.75  and  1 .5  (average 
TDT/Exp  =  1 .1  with  observed  incidence  of  8.9%).  Thus,  dives  with  ascent  rates  around 
30  fsw/min  have  incidence  almost  four  times  higher  than  dives  with  60  fsw/min  rates.  To 
test  the  hypothesis  that  slow  ascents  are  dangerous,  we  added  a  term  containing  the 
TDT/Exp  ratio  to  the  Equation  4  LOGIT  function.  The  result  indicates  that  a  30  fsw/min 
ascent  rate  may  actually  be  safer  than  a  60  fsw/min  rate,  but  the  result  is  not  statistically 
significant.  We  conclude  that  the  high  incidence  seen  near  a  ratio  of  1 .0  in  Figures  5D 
and  5E  probably  results  from  a  preponderance  of  aggressive  dive  profiles  (that  is,  dives 
that  are  long  and  deep)  in  the  30  fsw/min  group,  not  from  ascent  rate  itself. 

We  generate  Figure  6  from  the  No-Stop  calibration  data  by  forming  bins  according  to 
estimated  Pdcs  and  then  plotting  observed  DCS  for  a  bin  against  incidence  predicted  by 
the  P-No-Stop  Model  for  that  bin.  Solid  lines  and  curves  on  the  graphs  are  drawn  by 
least-squares;  equations  for  the  trends  and  R2  values  appear  in  boxes.  The  equations 
refer  to  probability  fractions  rather  than  to  percentages.  The  person-dives  are  clustered 
in  four  of  the  five  bins  with  DCS  incidence  less  than  5%,  and  three  of  the  bins  have 
observed  incidence  less  than  2%.  Observed  and  predicted  outcomes  are  highly 
correlated.  Chi-square  is  2.81  for  the  points  in  Figure  6A,  well  below  the  chance  value 
of  9.49  for  95%  confidence  for  4  degrees  of  freedom. 

In  the  low  probability  range  (Figure  6B),  agreement  is  fair:  for  the  data  point  with 
predicted  incidence  of  1 .26%,  observed  incidence  is  1 .23%,  but  the  other  points  deviate 
from  the  identity  line.  In  other  work  we  recommend  adjusting  or  fine-tuning  a 
probabilistic  model  when  the  low-incidence  range  appears  to  have  a  systematic  bias.24 
The  box  contains  the  equation  for  the  curve  drawn  by  least  squares  through  the  data 
points.  According  to  the  curve  and  to  this  equation,  observed  incidence  is  1.71%  when 
predicted  incidence  is  2.0%  and  observed  incidence  is  1.0  when  predicted  incidence  is 
0.64%.  The  differences  between  predicted  and  observed  incidences  in  Figure  6B 
appear  to  be  systematic,  so  we  judge  that  applying  the  fine-tuning  process  may  be 
appropriate  for  the  P-No-Stop  Model.  The  crux  of  the  process  is  to  substitute  the 
“probable  incidence,”  defined  as  the  average  observed  incidence  calculated  by  the 
least-squares  trend  curve,  for  the  predicted  incidence. 
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incidence  Pred 


Figure  6.  Observations  (DCS  cases/100  person-dives)  vs.  P-No-Stop  Model  predictions 
within  approximately  equal  subdivisions  of  the  calibration  data.  Dashed  lines  are  identity 
lines.  Horizontal  and  vertical  lines  show  95%  confidence  intervals  around  the  predicted 
and  observed  incidences.  A:  all  the  calibration  data;  the  solid  line  is  the  least-squares 
line  (see  the  equation  in  the  box).  B:  low-incidence  points  from  Panel  A;  the  solid  curve 
is  the  least-squares  curve  for  the  4  points  (see  the  equation  in  the  box). 


Goodness  of  fit,  subsaturation  dives 

Figure  7  shows  Pdcs  isopleths  in  the  subsaturation  range  (arbitrarily  defined  as  bottom 
times  between  240  and  1 ,440  min)  for  the  no-stop  probabilistic  model.  Clearly  the 
estimates  are  beyond  the  range  of  the  data  except  for  relatively  deep  subsaturation 
dives  and  for  dives  with  high  Pdcs.  There  are  no  dives  at  all  in  the  calibration  dataset 
for  depths  shallower  than  20  fswg,  and  no  DCS  cases  shallower  than  23  fswg.  Almost 
all  the  dives  in  the  box  have  high  estimated  Pdcs,  so  Pdcs  estimates  less  than  5%  tend 
to  be  extrapolations  from  regions  where  the  bulk  of  the  calibration  data  lie. 

The  two  dashed  vertical  line  segments  at  the  top  of  Figure  7  show  1  %  Pdcs  (left)  and 
2%  Pdcs  (right)  predicted  by  a  model  based  on  Equation  4  but  using  saturation  dives 
only.14  The  P-No-Stop  Model  curves  for  1  and  2%  Pdcs  match  up  well  with  their 
counterparts  from  the  other  model. 

For  an  example  of  the  relation  between  the  Pdcs  isopleths  and  DCS  incidences, 
consider  the  20%  Pdcs  curve  in  Figure  7.  The  curve  passes  through  two  groups  of 
triangles:  (1 )  a  group  with  bottom  times  of  1 ,440  min  with  depths  between  26  and  38 
fswg  represents  8  DCS  cases;  there  are  25  person-dives  at  1 ,440  min,  so  the  average 
incidence  is  32%,  and  (2)  a  group  with  720  min  bottom  time  between  31  and  41  fswg 
represents  23  DCS  cases;  there  are  118  person-dives  at  720  min,  so  average  incidence 
is  19.5%.  Thus,  the  groups  have  appropriate  incidences  to  be  near  the  20%  isopleth. 
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In  Figure  7  only  two  triangles  are  below  or  to 
the  left  of  the  isopleth  for  1  %  Pdcs.  One  is  at 
an  equivalent  depth  of  30  fswg  and  bottom 
time  of  240  min,  with  a  square  around  it  to 
signify  that  it  is  one  of  the  points  for 
breathing  gas  having  F02  above  that  of  air 
[case  number  (1)  in  Table  2].  The  other,  at 
the  lower-right  corner  of  the  graph,  is  for  two 
DCS  cases  superimposed  on  each  other  at 
an  equivalent  depth  of  57  fswg  with  bottom 
time  of  60  min.  The  circle  around  it  signifies 
that  it  is  for  breathing  of  F02  below  that  of  air 
[cases  (2)  and  (3)  in  Table  2],  Several 
triangles  are  below  or  to  the  left  of  the  curve 
for  2%  Pdcs. 

To  probe  the  fit  of  the  P-No-Stop  Model  in  the 
subsaturation  range,  we  sort  the  dives  in  the 
No-Stop  spreadsheet  by  depth  after 
excluding  dives  with  bottom  times  of  1 ,440 
min  or  longer  (saturation  dives).  From  the 
shallow  dives  in  the  sorted  set,  we  produce 
three  bins  centered  at  30,  35,  and  40  fswg.  The  bins  for  35  and  40  fswg  contain  enough 
person-dives  that  it  is  feasible  to  sort  them  again  by  bottom  time  to  make  additional 
bins.  The  Figure  8  list  of  averages  of  depths  and  bottom  times  for  the  seven  bins  that 
result  can  be  compared  with  Figure  7.  No  subsaturation  dives  are  shallower  than  30 
fswg  in  the  dataset.  The  Figure  8  graphs  show  that  the  No-Stop  probabilistic  model  fits 
well  with  the  limited  data  available:  the  chi-square  for  Figure  8A  is  1 .4,  and  12.6  is  the 
tabulated  value  for  95%  confidence  and  6  degrees  of  freedom. 

DETERMINISTIC  MODEL 

In  Figures  9,  10,  and  11 ,  the  traces  often  lie  close  to  each  other,  making  it  difficult  to 
distinguish  one  from  the  others.  This  is  an  important  point:  the  various  options  for  no¬ 
stop  dives  do  not  differ  from  each  other  significantly,  if  we  consider  the  inherent 
uncertainty  in  modeling  dive-outcome  data.  Figure  9  shows  how  the  models  differ  from 
each  other  in  an  important  way  only  with  deep,  short  dives.  The  top  panels  in  Figure  9 
show  symbols  for  dive-outcome  data  along  with  our  models’  options  for  no-stop 
instructions.  The  bottom  panels  illustrate  the  relationship  of  our  models  to  the  current 
U.S.  Navy  no-stop  instructions  and  the  uncertainty  involved  in  our  estimates. 

In  Figure  9Athe  beaded  trace  for  1%  Pdcs  according  to  our  No-Stop  probabilistic  model 
avoids  all  triangles  and  bull’s-eye  symbols  except  for  a  triangle  inside  a  large  square  at 
30  fswg.  As  in  Figure  3,  we  drew  the  solid  heavy  curve  for  our  deterministic  model  by 


Figure  7.  Probability  isopleths  for  the  P-No- 
Stop  Model  (solid  curves)  for  shallow  dives. 
The  dashed  box  outlines  the  region  of  sub¬ 
saturation  dives.  The  Pdcs  values  are,  from 
left  to  right,  1,2,3,  and  20%.  Symbols  are  as 
in  Figure  3.  Dashed  vertical  line  segments  at 
the  top  left  of  the  graph  are  for  1  and  2%  Pdcs 
according  to  a  model  based  only  on  saturation 
dives. 
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Figure  8.  Subsaturation  dives:  incidences  observed  and  predicted;  the  inset  lists  the 
divisions  of  the  data  plotted.  Sloping  solid  lines  are  trend  lines  with  accompanying 
equations.  Horizontal  and  vertical  lines  are  95%  confidence  intervals  around  the  predicted 
and  observed  DCS  incidences.  Dashed  lines  are  identity  lines.  A:  all  points  in  the  inset  list. 

B:  five  lowest  incidences  in  the  inset  list. 

eye;  it  tends  to  superimpose  on  the  1%  Pdcs  trace  for  the  P-No-Stop  Model  at  the  right 
of  Figures  9A  and  9C  and  also  at  the  left  of  Figures  9B  and  9D.  In  Figure  9B  the  middle 
part  of  the  D-No-Stop  Model  has  many  circles  below  it,  evidence  that  the  location  of  the 
D-No-Stop  curve  appropriately  divides  DCS  cases  from  non-DCS  cases. 

The  important  feature  of  Figure  9  is  that  the  deterministic  model  traces  deviate  from  the 
probabilistic  model  traces  at  the  right  of  Figures  9B  and  9D.  In  avoiding  triangles  and 
bull’s-eye  symbols  on  the  right  side  of  Figure  9B,  the  D-No-Stop  Model  allows 
considerably  less  bottom  time  than  the  No-Stop  curves  for  1%  and  2%  Pdcs.  The  DCS 
cases  that  cast  doubt  on  the  safety  of  the  beaded  P-No-Stop  Model  prescriptions  at  the 
right-hand  side  of  Figure  9B  are  the  lowest  bull’s-eye  point  (165  fswg,  bottom  time  of  10 
min,  2  DCS  cases)  and  two  triangles  (198  fswg,  bottom  time  of  9  min,  1  DCS  case;  and 
231  fswg,  bottom  time  of  6.8  min,  2  DCS  cases). 
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Figure  9.  The  deterministic  D-No-Stop  Model  (heavy  traces  that  go  no  higher  than  140  min) 
compared  with  isopleths  for  1  and  2%  Pdcs  (beaded  traces)  for  the  probabilistic  P-No-Stop 
Model.  Symbols  are  as  in  Figure  3.  A  and  C:  shallow  dives;  B  and  D:  deep  dives.  A:  the 
dashed  box  shows  the  bottom  of  the  subsaturation  range  seen  in  Figure  7;  at  the  lower  right, 
the  probabilistic  No-Stop  Model’s  1%  Pdcs  curve  is  obscured  by  the  D-No-Stop  curve.  B:  for 
deep  dives,  the  deterministic  model  deviates  from  two  traces  for  the  probabilistic  model.  C: 
the  95%  confidence  intervals  (gray)  around  the  probabilistic  model’s  estimate  of  1%  Pdcs; 
the  trace  with  black  squares  is  for  USN57;  the  probabilistic  model’s  trace  for  2%  Pdcs  is 
omitted  to  reduce  clutter.  D:  as  in  panel  C,  for  deep  dives. 

The  bottom  times  permitted  by  the  P-No-Stop  Model  for  depths  between  140  and  190 
fswg  are  as  much  as  twice  as  long  as  those  for  the  D-No-Stop  Model.  The  USN57 
trace,  also  deterministic,  is  close  to  the  No-Stop  2%  curve  at  the  left  of  Figure  9D  but 
falls  below  the  1  %  trace  at  the  right. 

The  1%  curve  for  the  P-No-Stop  Model  and  the  D-No-Stop  Model  give  similar 
prescriptions  where  the  probabilistic  model  is  successful  in  Figure  9,  at  depths  from  40 
through  about  130  fswg.  However,  for  no-stop  dives  deeper  than  1 30  fswg  and  1% 
Pdcs  (right  side  of  Figures  9B  and  9D),  the  deterministic  model  calls  for  much  shorter 
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bottom  times  than  the  P-No-Stop  Model.  Our  interpretation  is  that  the  probabilistic 
model  fails  for  deep  no-stop  dives  whereas  the  deterministic  model,  which  avoids  DCS 
cases,  is  appropriate.  Furthermore,  we  speculate  that  Pdcs  for  the  D-No-Stop  Model 
may  be  about  1  %  for  dives  deeper  than  130  fswg. 

The  uncertainty  in  our  modeling  is  illustrated  by  the  95%  confidence  intervals  in  Figures 
9C  and  9D,  which  cover  a  span  that  is  about  the  same  as  the  difference  between  the 
different  curves  in  Figures  9A  and  9B.  In  Figure  9D  the  D-No-Stop  trace  just  crosses 
the  lower  confidence  boundary  at  the  extreme  right.  However,  the  confidence  interval 
boundaries  do  not  include  the  USN57  trace  at  the  right  of  Figure  9D. 

We  generated  our  model  for  standard  air  diving,  the  StandAir  Model,13  from  a  large 
calibration  dataset:  45%  of  the  dives  were  no-stop  dives,  essentially  those  in  the 
calibration  dataset  for  our  present  models,  and  55%  were  dives  with  decompression 
stops.  Figure  1 0  compares  traces  in  Figure  9  with  traces  for  the  StandAir  Model  and  for 
a  deterministic  model  known  as  the  VVal-1 8  Algorithm,25,26  which  gives  TDTs  for  dives 
having  decompression  stops  similar  to  those  for  our  StandAir  Model. 


Figure  10.  Traces  for  the  StandAir  Model  and  the  VVal-1 8  Algorithm  added  to  those  seen  in 
Figure  9.  A:  shallow  dives;  deeper  than  40  fswg,  the  probabilistic  P-No-Stop  Model’s  1%  Pdcs 
curve  is  partly  obscured  by  the  D-No-Stop  curve.  B:  deep  dives. 

•  In  Figure  1 0A  the  2%  Pdcs  trace  for  the  StandAir  Model  and  the  VVal-1 8  trace  lie 
between  the  1%  Pdcs  trace  (hidden  at  the  right  by  the  D-No-Stop  trace)  and  2% 
Pdcs  traces  for  the  P-No-Stop  Model  (beaded)  and  StandAir  Model  (diamond 
symbols). 

•  At  the  left  of  Figure  1 0B,  the  gray  trace  for  the  VVal-1 8  Algorithm  corresponds  to 
about  2%  Pdcs,  according  to  the  P-No-Stop  Model,  and  the  StandAir  Model  falls  to 
the  left  of  the  P-No-Stop  Model  for  1  %  Pdcs.  At  the  right  of  Figure  1 0  B,  the  VVal-1 8 
trace  is  slightly  below  the  probabilistic  traces,  and  the  trace  for  the  USN57 
instructions  is  far  below  all  the  other  traces. 
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•  The  most  important  aspect  of  Figure  10B  is  at  the  right-hand  side:  the  2%  trace  for 
the  Stand  Air  probabilistic  model  falls  between  the  two  beaded  traces  for  our  No-Stop 
probabilistic  model  and  well  above  the  trace  for  our  D-No-Stop  Model. 


DISCUSSION 

NO-STOP  LIMITS  OF  OTHER  NAVIES 

Table  4  lists  current  no-stop  instructions  for  operational  tables  used  by  navies  of  Great 
Britain,10  Canada,11  France,12  and  the  United  States.1  The  USN57  tends  to  allow  longer 
bottom  times  than  the  other  navies’  tables,  which  tend  to  be  similar  to  each  other. 
Information  from  our  two  models  appears  at  the  right  of  Table  4. 


U  =  unlimited  time  at  depth  is  allowed 
-  =  no  entry 


Using  the  format  introduced  in  previous  figures  (shallow  no-stop  dives  in  left-hand 
panels  and  deep  no-stop  dives  in  right-hand  panels),  Figure  11  compares  the 
information  in  Table  4.  In  Figure  11  A,  the  USN57  trace  is  a  little  below  the  2%  Pdcs 
curve  for  the  P-No-Stop  Model.  Except  for  USN57,  the  navy  traces  tend  to  be  a  little 
above  or  on  the  1%  Pdcs  trace;  apparently,  the  no-stop  prescriptions  of  the  British, 
Canadian,  and  French  navies  reflect  risks  near  1%  for  depths  shallower  than  100  fswg. 
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Figure  11.  No-stop  times  for  various  navies  plus  our  deterministic  D-No-Stop  Model  (heavy 
trace)  and  isopleths  for  our  No-Stop  probabilistic  model  (beaded  curves  for  2%  Pdcs  (upper) 
and  1%  Pdcs  (lower)).  In  both  panels,  USN57  is  the  trace  with  black  squares.  A:  shallow  dives. 
B:  deep  dives. 


For  deep  dives,  the  differences  between  our  P-No-Stop  Model  and  the  navy  instructions 
in  Figure  11 B  are  large;  all  the  navy  traces  mandate  much  shorter  bottom  times  than  do 
the  beaded  traces  for  our  probabilistic  model.  The  D-No-Stop  trace  is  between  the 
traces  for  the  other  navies  and  the  traces  for  the  probabilistic  model.  For  example, 

Table  4  and  the  right  side  of  Figure  11 B  show  that  the  allowed  bottom  times  for  no-stop 
dives  to  190  fswg  are  10  and  12  min  for  the  probabilistic  P-No-Stop  Model,  7.5  min  for 
the  D-No-Stop  Model,  and  5  min  for  USN57  and  Canada;  dives  to  190  fswg  are  not 
allowed  by  British  and  French  navies. 

Figure  12  shows  that  estimates  of  Pdcs  calculated  with  our  P-No-Stop  Model  tend  to  be 
1 .5%  to  2%  for  the  USN57  schedule  for  depths  shallower  than  100  fswg.  The 
uncertainty  of  the  USN57  estimates  is  illustrated  by  the  two  varying  dashed  traces  for 
upper  and  lower  95%  confidence  intervals.  Another  aspect  of  uncertainty  involves  the 
two  dashed  horizontal  traces:  the  least-squares  trend  curve  in  Figure  6  shows  that 
when  Pdcs  is  1%,  average  observed  incidence  or  “probable  incidence”  is  0.64%,  and 
when  Pdcs  is  2%,  probable  incidence  is  1.71%.  Thus,  the  probable  incidence 
associated  with  the  upper  horizontal  line  in  Figure  12  is  1 .71%,  and  the  probable 
incidence  for  the  lower  horizontal  trace  is  0.64%. 

In  Figure  12  the  other  navies  tend  to  have  Pdcs  around  1%,  or  less  than  2%  for  depths 
shallower  than  100  fswg,  and  confidence  intervals  around  their  traces  (not  shown)  are 
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Figure  12.  Pdcs  estimates  for  no-stop  times, 
calculated  by  our  No-Stop  probabilistic  model. 
Irregular  dashed  traces  show  the  95%  confidence 
intervals  for  USN57.  The  dashed  horizontal  line 
segments  show  “probable  incidences,"  which  are 
less  than  the  Pdcs  for  the  same  location. 


similar  to  the  intervals  for  the  USN57 
Model.  If  the  P-No-Stop  Model  is  faulty 
for  deep  dives,  its  estimates  of  Pdcs  for 
the  right  side  of  Figure  1 2B  are  incorrect. 


PITFALLS  IN  PROBABILISTIC 
MODELING 

Figures  9B,  9D,  10B,  and  1 1 B  illustrate  a 
discrepancy  between  deterministic  and 
probabilistic  models:  the  bottom  times 
prescribed  for  deep,  no-stop  dives  by  the 
deterministic  D-No-Stop  Model  and  the 
operational  navies  are  considerably 
shorter  than  the  bottom  times  prescribed 
by  the  probabilistic  P-No-Stop  and 
StandAir  models.  Our  interpretation  is 
that  the  probabilistic  models  are  at  fault. 
We  recognize  three  general  pitfalls  in 
probabilistic  modeling. 


Inappropriate  equations 

The  first  pitfall  is  that  the  equation  used  to  fit  the  data  may  introduce  bias.  The  simple 
analogy  here  is  that  it  is  undesirable  to  fit  a  least-squares  straight  line  to  points  that 
appear  as  a  curve  on  a  graph,  and  any  applications  of  the  parameter  values  derived 
from  such  a  fit  are  misleading.  The  inappropriate  equation  prevents  the  model  from 
providing  a  good  fit  to  the  data  points  in  certain  regions. 


Our  use  of  Equation  4  presumes  that  the  combination  of  two  compartments,  one  of 
them  weighted,  acceptably  represents  a  continuum  of  many  critical  compartments 
having  varying  halftimes.  However,  the  halftimes  of  21 .5  and  420  min  found  by  the 
statistical  process  can  be  expected  to  provide  questionable  results  for  deep  dives;  for 
such  dives,  bottom  times  are  comparatively  short,  so  the  model  should  probably  include 
a  compartment  with  halftime  shorter  than  21  min.  When  we  generated  a  three- 
exponential  model  with  two  of  the  exponential  terms  weighted,  one  of  the  halftimes  was 
11  min.  We  also  tested  models  with  a  third  exponential  term  having  halftimes  fixed  at  6, 
3,  or  1  min.  However,  for  all  these,  the  LL  was  not  better  than  that  for  Equation  4,  and 
the  prescriptions  for  deep,  short  dives  were  not  significantly  better  than  those  of  the  P- 
No-Stop  Model.  Apparently  the  lack  of  deep,  short  dives  in  the  dataset  is  more 
troublesome  than  is  the  lack  of  a  compartment  with  a  short  halftime. 
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Out-of-bounds  data 


A  second  pitfall  in  probabilistic  modeling  concerns  the  quality  of  the  calibration  data. 

The  practical  value  of  the  model  is  compromised  if  the  data  are  not  representative  of  the 
intended  use  of  the  model.  Our  calibration  dataset  provides  an  unavoidable  example  of 
this  defect.  The  P-No-Stop  Model  is  intended  for  use  with  dives  having  the  current  U.S. 
Navy  ascent  rate  of  30  fsw/min,  but  about  half  the  calibration  dives  have  an  ascent  rate 
of  60  fsw/min;  these  60  fswg/min  dives  are  “out-of-bounds”  data.  If  rapid  ascents  are 
more  hazardous  than  slow  ascents,  the  resulting  bias  is  on  the  side  of  safety.  Some  of 
the  saturation  dives  have  ascent  rates  slower  than  30  fsw/min,  but  this  slow  rate  may 
not  be  troublesome  if  the  DCS  bubbles  that  occur  after  saturation  dives  are  only  in 
compartments  with  very  slow  halftimes. 

Including  data  beyond  the  range  of  the  desired  outcome  invites  error:  the  out-of-bounds 
data  may  distort  the  fit  in  the  desired  region,  in  this  report  our  region  of  interest  includes 
no-stop  dives  ranging  from  the  few  minutes  of  bottom  time  that  are  appropriate  for  a 
1 90-fswg,  no-stop  dive  to  the  days-long  bottom  times  for  saturation  dives.  Such  a  wide 
region  of  interest  has  merit  for  deriving  prescriptions  for  subsaturation  dives,  where  data 
are  scarce,  because  the  subsaturation  region  is  interpolated  between  standard  air  dives 
and  saturation  dives. 

In  another  study  we  tested  our  contention  that  the  calibration  data  used  to  produce  a 
model  should  be  limited  to  the  types  of  dives  with  which  the  model  is  to  be  used.27  We 
produced  a  model  that  included  both  saturation  and  nonsaturation  dives  with,  in  both 
cases,  decompression  stops.  For  prescriptions  of  TDTs,  this  combination  model  was 
clearly  inferior  to  an  alternative  model  based  on  a  dataset  that  excluded  saturation 
dives. 

This  report  allows  a  second  test  of  our  contention  that  out-of-bounds  data  can  bias  a 
model.  Figure  10  shows  curves  for  1%  and  2%  Pdcs  for  no-stop  dives  according  to  our 
previous  StandAir  Model,13  which  was  generated  from  a  database  that  contained 
approximately  equal  numbers  of  no-stop  dives  and  dives  with  decompression  stops. 
Thus,  so  far  as  no-stop  dives  are  concerned,  the  StandAir  Model  contains  out-of-bounds 
dives  having  decompression  stops.  However,  the  StandAir  Model  prescribes  essentially 
the  same  no-stop  depths  and  bottom  times  for  2%  Pdcs  as  our  P-No-Stop  Model 
prescribes.  This  agreement  between  models  indicates  that  including  the  out-of-bounds 
dives  that  have  decompression  stops  does  not  unduly  bias  the  no-stop  predictions  of 
the  StandAir  Model.  The  agreement  between  the  two  models  is  not  as  good  for  Pdcs  of 
1  %  as  it  is  for  Pdcs  of  2%. 


Balance 

The  third  pitfall  concerns  the  balance  of  the  calibration  data.  The  poor  balance  of  the 
available  data  is  illustrated  in  Figures  2  and  5.  In  Figure  2A,  person-dives  are  scarce, 
and  DCS  incidence  is  low  for  deep  dives.  Figure  5A,  with  bins  containing  approximately 
equal  numbers  of  person-dives,  gives  another  perspective:  the  incidence  of  DCS  is  low 
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for  dives  around  50  to  100  fswg.  The  poor  distribution  of  bottom  times  is  manifested  in 
Figures  2B,  5B,  and  5C.  Figure  2B  shows  that  the  numbers  of  person-dives  (main 
figure)  and  incidence  (inset)  differ  widely  between  bins;  Figure  5B  has  a  peak  of  DCS 
incidence  for  saturation  dives;  and  Figure  5C  has  another  peak  around  40  min  and  very 
low  incidence  below  30  min.  Figure  2C  shows  that  ascent  rate,  characterized  by  the 
TDT/Exp  ratio,  is  very  unevenly  distributed,  especially  with  respect  to  DCS  incidence 
(inset  of  Figure  2C  and  Figures  5D  and  6E),  but  ascent  rate  does  not  directly  affect  the 
P-No-Stop  Model  because  Equation  4  does  not  include  ascent  rate  or  TDT  as  a 
variable.  Although  thousands  of  dives  are  recorded  in  the  U.S.  Navy  Decompression 
Database,4,5  our  models  are  weakened  by  lack  of  data  in  particular  depth/bottom-time 
regions  such  as  saturation  dives  (only  246  dives  with  22  DCS  cases,  none  at  depths 
shallower  than  23  fswg),  subsaturation  dives  (no  dives  at  all  between  720  and  1 ,400 
min  bottom  times),  and  deep  dives. 

Data  determine  the  outcome  of  the  statistical  process,  of  course,  and  depth/bottom-time 
regions  that  include  many  data  points  have  more  impact  on  the  model’s  parameters, 
and  therefore  on  the  model’s  prescriptions,  than  do  regions  that  include  few  points.  It 
follows  that  the  probabilistic  fitting  of  poorly  balanced  data  may  give  faulty  estimates  in 
some  regions  of  interest.  In  regions  where  data  are  scarce,  the  model’s  output  is  biased 
by  regions  that  contain  more  data.  If  no  data  are  in  a  region,  the  model’s  output  for  that 
region  is  interpolated  or  extrapolated  from  regions  where  dive  outcome  information  is 
plentiful.  Outer  ranges  of  a  variable  —  such  as  the  shallowest  and  deepest  dives,  dives 
with  the  shortest  or  longest  bottom  times,  and  dives  with  the  lowest  or  highest  Pdcs  — 
are  especially  vulnerable  to  suboptimal  fitting  equations  and  poorly  balanced  data. 
Accordingly,  our  results  show  that  the  Pdcs  estimates  given  by  our  P-No-Stop  Model  are 
questionable  for  the  deepest  dives. 

A  partial  remedy  for  lack  of  data  is  to  seek  guidance  from  outside  sources.  One  type  of 
outside  source  is  the  operational  navy  tables  discussed  under  the  heading  NO-STOP 
LIMITS  OF  OTHER  NAVIES;  the  navies  prescribe  shorter  bottom  times  for  deep  dives 
than  our  probabilistic  model  prescribes.  A  second  outside  source  of  information  is  our 
use  of  selected  individual  dives  with  decompression  stops  (bull’s-eyes  in  the  figures) 
with  our  deterministic  models  to  attempt  a  partial  compensation  for  the  paucity  of  DCS 
cases  in  deep  no-stop  dives. 

TABLES  OF  DCS  PROBABILITY 

In  Tables  5  through  8  our  recommendations  appear  in  heavy  boxes.  Information  outside 
the  heavy  boxes  allows  comparisons  between  models.  The  recommendations  are 
essentially  arbitrary  and  based  on  presently  available  data;  future  developments  may 
modify  or  reverse  them. 

■  We  recommend  prescriptions  and  Pdcs  estimates  given  by  the  P-No-Stop  Model 
for  subsaturation  dives  and  for  standard  dives  shallower  than  1 40  fswg. 
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■  For  no-stop  dives  deeper  than  1 30  fswg,  we  recommend  the  prescriptions  of  the 
D-No-Stop  Model  —  which,  we  speculate,  represents  Pdcs  near  1  %.  We  judge 
that  our  probabilistic  model  does  not  give  acceptable  prescriptions  and 
appropriate  Pdcs  estimates  deeper  than  1 30  fswg. 

Standard  air  no-stop  dives:  Table  5 

The  upper  heavy  box  in  Table  5  lists  no-stop  bottom  times  for  Pdcs  of  1  %,  2%,  and  3% 
with  an  ascent  rate  of  30  fswg/min,  according  to  our  P-No-Stop  Model.  The  2%  Pdcs 
no-stop  times  in  Table  5  are  longer  than  the  USN57  no-stop  bottom-time  limits  in  force 
in  the  year  2004  for  all  depths  except  35  fswg.  The  USN57  bottom  times  are  in  the 
range  of  1%  and  2%  Pdcs,  according  to  the  P-No-Stop  Model,  except  for  dives  to  110 
fswg  and  deeper,  where  the  USN57  times  are  shorter  than  those  given  for  1%  Pdcs  in 
the  table. 


TABLE  5.  NO-STOP  TIME  LIMITS  (MIN):  THE  U.S.  NAVY  DIVING 
TABLE  (USN57),  THE  P-NO-STOP  MODEL,  AND  THE  D- 
NO-STOP  MODEL 

Heavy  boxes  indicate  our  recommended  bottom  times. 

Light  solid  boxes  =  bottom  times  less  than  those  for  USN57. 

Light  print  =  questionable  values. 

(1)  (2)  (3)  (4)  (5)  (6)  (7) 

Depth, 

fswg 

USN57 

<-■  No-Stop  Model  -> 

1%  2%  3% 

D-No-Stoo 

Pdcs  for 
D-No-Stoo 

25 

595 

462 

694 

889 

450 

1. 0% 

30 

405 

294 

429 

527 

300 

1.0% 

35 

310 

198 

293 

358 

200 

1.0% 

40 

200 

137 

208 

257 

137 

1.0% 

50 

100 

75 

112 

141 

76 

1 .0% 

60 

60 

51 

70 

86 

51 

1.0% 

70 

50 

38 

51 

60 

38.5 

1.0% 

80 

40 

31 

40 

46 

31 

1.0% 

90 

30 

26 

33 

37 

26 

1.0% 

100 

25 

22 

28 

32 

22 

0.9% 

110 

20 

20 

24 

27 

19 

0.9% 

120 

15 

18 

22 

24 

17 

0.9% 

130 

10 

16 

19 

22 

15 

0.8% 

140 

10 

15 

18 

20 

13  | 

150 

5 

13 

16 

18 

11.5 

— 

160 

5 

12 

15 

16 

10 

— - 

170 

5 

12 

14 

15 

9 

— 

180 

5 

11 

13 

14 

8 

— 

190 

5 

10 

12 

13 

7.5 

28 


For  dives  deeper  than  130  fswg,  we  believe  our  No-Stop  probabilistic  model  is  too 
liberal;  the  bottom  times  in  the  heavy  box  in  Table  5,  column  6,  are  from  the  D-No-Stop 
Model.  Even  these  are  as  much  as  two  and  a  half  times  longer  than  times  for  USN57. 
The  largest  difference  is  at  150  fswg,  where  the  time  is  11 .5  min  instead  of  the  5  min  of 
USN57.  We  recommend  caution  in  lengthening  no-stop  bottom  times  for  deep  dives, 
because  such  dives,  at  least  in  animal  studies,  have  been  associated  with  DCS  signs 
and  symptoms  that  are  more  serious  than  those  occurring  in  shallower  dives.28-30 
However,  when  we  scrutinize  information  about  individual  DCS  cases  in  the  No-Stop 
spreadsheet,  we  cannot  confirm  a  trend  for  deep  cases  to  be  more  serious  than  shallow 
cases. 

If  we  read  horizontally  along  bottom-time  entries  for  dives  deeper  than  1 30  fswg  in 
columns  3  to  5  of  Table  5,  we  see  that  one  to  three  minutes  of  change  in  bottom  time 
causes  a  1%  change  in  Pdcs.  Therefore,  if  the  entries  in  column  6  represent  1%  Pdcs, 
entries  for  dives  deeper  than  140  fswg  in  column  4  may  correspond  to  Pdcs  as  high  as 
5%. 

Standard  air  dives  over  a  range  of  risk:  Table  6 

Table  6,  which  is  similar  to  a  Table  1  published  in  the  description  of  the  NMRI  ‘93 
probabilistic  model,15  shows  how  bottom  times  increase  as  Pdcs  increases.  Two 
columns  near  the  center  of  Table  6  show  that  the  P-No-Stop  Model’s  bottom  times  for 
2.3%  Pdcs  are  longer  than  NMRI  ‘93’s  for  the  same  Pdcs,  more  than  50%  longer  for 
shallow  dives.  At  the  bottom  of  Table  6,  the  second  set  of  entries  for  depths  from  1 40 
through  190  fswg  show  the  bottom  times  obtained  from  the  D-No-Stop  model. 

Differences  between  the  basic  properties  of  the  NMRI  ’93  Model  and  our  probabilistic 
model  are  revealed  by  dose-response  curves.  In  Figure  13,  traces  for  our  P-No-Stop 
Model  have  sharper  elbows  than  traces  for  the  NMRI  ‘93  Model,  so  that  the  curves  for 
the  two  models  tend  to  cross  each  other.  The  horizontal  line  segments  in  each  panel 
show  the  2.3%  level  of  risk  deemed  acceptable  by  developers  of  the  NMRI  ‘93  Model. 
The  horizontal  line  segment  in  the  40  fswg  panel  shows  that  the  P-No-Stop  Model 
allows  a  considerably  longer  bottom  time  than  the  NMRI  ‘93  Model  allows  for  the  same 
Pdcs,  and  the  vertical  line  segment  shows  that  NMRI  ’93  requires  a  bottom  time 
equivalent  to  only  about  1%  risk,  according  to  the  P-No-Stop  Model.  The  crossovers  in 
the  other  panels  cause  the  NMRI  ‘93  Model  to  call  for  shorter  bottom  times  than  the  P- 
No-Stop  Model  at  low  Pdcs,  but  for  longer  bottom  times  at  high  Pdcs. 

Except  at  40  fswg,  the  traces  for  the  two  models  tend  to  be  close  to  each  other  when 
Pdcs  is  2  to  4%  in  Figure  1 3.  The  vertical  line  in  the  1 80-fswg  panel  is  for  8  min,  the 
time  our  D-No-Stop  Model  prescribes  for  1 80  fswg.  The  NMRI  ‘93  Model  curve  meets 
the  heavy  vertical  line  at  about  2%  Pdcs,  but  the  P-No-Stop  Model  meets  the  vertical 
line  at  a  very  low  Pdcs.  If  the  D-No-Stop  Model’s  vertical  line  segment  is  the  correct 
bottom  time  for  1%  Pdcs  at  180  fswg,  the  NMRI  ’93  Model  overestimates  the  risk  and 
our  P-No-Stop  Model  underestimates  it. 


TABLE  6.  BOTTOM  TIMES  (MIN)  ACCORDING  TO  THE  P-NO-STOP  MODEL 
FOR  VARIOUS  DCS  PROBABILITIES 

Boldface  columns  show  current  U.S.  Navy  no-stop  limits  (left)  and  the  recommendations  of 
the  NMRI  ‘93  model  (middle). 

Light  print  =  questionable  values . 

For  depths  from  140  through  190  fswg,  we  recommend  the  D-No-Stop  Model  bottom  times 
listed  in  the  heavy  box  at  the  bottom  of  the  table. 


USN57  1.0%  2.0%  2.1%  2.2%  2.3% 

Depth 


120 

15  1 

130 

10  I 

140 

i°  T 

150 

5 

160 

5 

170 

5 

180 

5 

190 

5 

140 

10 

150 

5 

160 

5 

170 

5 

180 

190 

5 

NMRI 

‘93  2.4%  2.5%  3.0%  4.0%  5.0%  6.0% 

2.3% 


774  793  889 

471  481  527 

321  328  358 
229  234  257 
167  171  188 

125  127  141 
96  98  108 

77  78  86 

55  56  60 

42  43  46 

35  35  37 

30  30  32 

26  26  27 

23  23  24 

20  21  22 


**  =  more  than  2,880  min 
-  =  no  entry 
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Bottom  time,  min  Bottom  time,  min 


Figure  13.  Dose-response  curves  for  the  P-No-Stop  Model 
(solid  curves)  and  the  NMRI  ‘93  Model  (dashed  curves)  for  four 
depths.  Horizontal  line  segments  show  the  no-stop  times  for 
NMRI  ‘93.  The  vertical  dashed  line  segment  in  the  40-fswg 
panel  compares  the  two  models  at  the  same  bottom  time.  The 
vertical  line  segment  in  the  180-fswg  panel  shows  the  D-No- 
Stop  Model’s  bottom  time. 


Subsaturation  dives:  Tables  7 
and  8 


Table  7  lists  bottom  times  for 
subsaturation  dives  calculated 
by  the  P-No-Stop  model  for  1 
and  2%  Pdcs.  Practical 
experience  indicates  that  the 
USN57  table  is  reliable  in  the 
subsaturation  region.31  For 
2%  Pdcs,  the  P-No-Stop 
Model’s  bottom  times  are  1 7% 
longer  than  those  of  USN57  at 
25  fswg,  5%  longer  at  30  fswg, 
and  6%  shorter  at  35  fswg. 

The  upper  and  lower  95% 
confidence  intervals  in 
columns  4,  5,  7,  and  8  give  an 
idea  of  the  uncertainty  involved 
in  the  bottom  times:  the 
confidence  intervals  extend 
over  100  minutes  on  either 
side  of  the  prescribed  times  at 
the  top  of  the  table;  for  40  fswg 
at  the  bottom  of  the  table, 
confidence  intervals  extend 
about  a  half  hour  from  either 
side  of  the  prescribed  times. 


Table  8  lists  bottom  times  less  than  2  days  given  by  the  P-No-Stop  Model  for 
subsaturation  and  shallow  dives  between  18  and  60  fswg  for  various  Pdcs  values.  For 
example,  the  2%  Pdcs  column  shows  that  a  bottom  time  of  1 ,735  min  (29  hr)  at  20  fswg 
and  a  bottom  time  of  694  min  (11 .6  hr)  at  25  fswg  have  the  same  risk  of  DCS.  The 
USN57  column  allows  405  min  (less  than  7  hr)  at  30  fswg,  a  bottom  time  that  gives  a 
little  less  than  2%  Pdcs.  The  1%  column  shows  that  a  diver  can  remain  at  35  fswg  for 
197  min  (3.3  hr),  and  the  5%  column  shows  that  a  diver  can  remain  at  35  fswg  for  453 
min  (7.6  hr).  At  35  fswg,  Pdcs  is  greater  than  30%  for  a  diver  who  makes  a  no-stop 
ascent  after  1 ,440  min  (24  hr). 
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TABLE  7.  BOTTOM  TIMES  (MIN)  FOR  SUBSATURATION  NO- 

STOP 

DIVES 

Light 

print  = 

interpolations  between  standard  air  dives  and 

!  saturation  dives  due  to  lack  of  data. 

Entries  in  light  boxes  =  bottom  times  less  than  those  for  USN57. 

(1) 

(2) 

(3) 

(4)  . 

(5) 

(6) 

(7) 

(8)  . 

No-Stop 

No-Stop 

Depth, 

Model, 

Low 

High 

Model, 

Low 

High 

fswg 

USN57 

1  %  Pdcs 

C.I. 

C.l. 

2%  Pdcs 

C.l. 

C.l. 

23 

572 

433 

693 

897 

715 

1,083 

24 

512 

391 

616 

783 

630 

928 

25 

595 

462 

353 

552 

694 

566 

811 

26 

419 

322 

499 

622 

508 

722 

27 

381 

292 

452 

562 

462 

648 

28 

348 

266 

411 

511 

422 

587 

29 

319 

244 

377 

468 

387 

535 

30 

405 

294 

222 

347 

429 

357 

490 

31 

270 

205 

320 

396 

329 

451 

32 

249 

188 

297 

366 

304 

416 

33 

230 

172 

275 

339 

282 

386 

34 

213 

159 

255 

315 

260 

359 

35 

310 

197 

145 

238 

293 

242 

334 

40 

200 

137 

98 

168 

208 

167 

241 

CONCLUSIONS 

1 .  Our  No-Stop  probabilistic  model,  based  on  all  the  no-stop  data,  gives  acceptable 
Pdcs  estimates  for  no-stop  dives  shallower  than  130  fswg,  but  not  for  no-stop 
dives  deeper  than  1 30  fswg. 

2.  For  deep  no-stop  air  dives,  our  D-No-Stop  model,  fashioned  to  avoid  observed 
cases  of  DCS,  gives  acceptable  bottom  times:  Pdcs  is  indeterminate,  but  we 
judge  it  to  be  about  1%. 

3.  The  P-No-Stop  Model  seems  to  give  acceptable  bottom-time  limits  for 
subsaturation  dives,  although  the  model  results  are  interpolations  between 
standard  air  dives  and  saturation  dives  in  much  of  this  region. 

4.  Our  evidence  that  the  P-No-Stop  Model  underestimates  risk  of  DCS  in  deep 
dives  shows  that  a  probabilistic  model  that  serves  well  for  much  of  the 
depth/bottom-time  region  of  the  dataset  can  give  faulty  results  in  a  certain  region, 
even  when  such  a  model  is  based  on  a  relatively  homogeneous  dataset 
containing  only  no-stop  dives. 
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TABLE  8.  SUBSATURATION  DIVES  AND  SHALLOW  DIVES:  ESTIMATES  OF 
BOTTOM-TIME  LIMITS  (MIN)  FROM  THE  P-NO-STOP  MODEL 

+++  =  bottom  times  of  48  hr  or  more 

Entries  in  light  boxes  =  bottom  times  less  than  those  for  USN57 

Light  print  =  interpolations  between  standard  air  dives  and  saturation  dives. 


—  Bottom  times,  min  — 
3%  5%  10% 


30%  50% 


1,293 

1,026 

858 

737 

646 

572 

512 

462 

419 

381 

348 

319 

294 

270 

249 

230 

213 

197 

183 

170 

158 

147 

137 

127 

119 

112 

105 

99 

93 

88 

83 

79 

75 

72 

69 

66 

63 

61 

58 

56 

54 

53 

51 


+++ 
+++ 
1,735 
1,288 
1,053 
897 
783 
694 
622 
562 
511 
468 
429 
396 
366 
339 
315 
293 
273 
255 
238 
222 
208 
195 
183 
171 
161 
151 
142 
134 
126 
119 
112 
106 
101 
96 
91 
87 
83 
80 
76 
73 
I  70 


+++ 

+++ 

+++ 

2,556 

1,561 

1,227 

1,028 

889 

784 

701 

633 

576 

527 

484 

447 

414 

385 

358 

334 

312 

292 

274 

257 

241 

226 

213 

200 

188 

178 

167 

158 

149 

141 

133 

126 

120 

114 

108 

103 

98 

94 

89 

86 


+++ 

+++ 

+++ 

+++ 

+++ 

2,798 

1,644 

1,294 

1,088 

945 

836 

750 

680 

621 

570 

526 

487 

453 

422 

394 

369 

346 

325 

306 

288 

271 

256 

242 

228 

216 

204 

193 

183 

173 

164 

156 

148 

141 

134 

127 

121 

115 

110 


+++ 

+++ 

+++ 

+++ 

+++ 

+++ 

+++ 

+++ 

2,577 

1,664 

1,334 

1,134 


993 

1,799 

+++ 

+++ 

886 

1,435 

+++ 

+++ 

800 

1,221 

1,888 

+++ 

729 

1,072 

1,499 

+++ 

669 

958 

1,275 

+++ 

617 

868 

1,120 

2,135 

572 

794 

1,004 

1,642 

533 

731 

911 

1,386 

497 

677 

834 

1,214 

466 

629 

769 

1,087 

437 

.  588 

714 

987 

411 

550 

665 

905 

387 

517 

622 

836 

365 

487 

584 

777 

345 

459 

549 

725 

326 

434 

518 

680 

309 

411 

490 

639 

293 

390 

464 

603 

278 

370 

440 

570 

264 

352 

418 

540 

251 

334 

398 

512 

238 

318 

379 

487 

227 

303 

361 

464 

215 

289 

344 

442 

205 

276 

329 

422 

195 

264 

314 

403 

186 

252 

300 

385 

177 

241 

287 

369 

169 

230 

275 

353 

161 

220 

263 

339 

154 

211 

252 

325 

33 


PAGE  INTENTIONALLY  LEFT  BLANK 


REFERENCES 


1.  Commander,  Naval  Sea  Systems  Command,  U.S.  Navy  Diving  Manual,  Revision 
4,  Publication  SS521-AG-PRO-010  (Arlington,  VA:  Naval  Sea  Systems 
Command,  1999). 

2.  P.  K.  Weathersby,  L.  D.  Homer,  and  E.  T.  Flynn,  “On  the  Likelihood  of 
Decompression  Sickness,”  JAppI  Physiol,  Vol.  57  (1984),  pp.  815-824. 

3.  E.  C.  Parker,  S.  S.  Survanshi,  P.  B.  Massell,  and  P.  K.  Weathersby,  “Probabilistic 
Models  of  the  Role  of  Oxygen  in  Human  Decompression  Sickness,”  JAppI 
Physiol,  Vol.  84  (March  1998),  pp.  1096-1102. 

4.  D.  J.  Temple,  R.  Ball,  P.  K.  Weathersby,  E.  C.  Parker,  and  S.  S.  Survanshi,  The 
Dive  Profiles  and  Manifestations  of  Decompression  Sickness  Cases  after  Air  and 
Nitrogen-Oxygen  Dives,  NMRC  Technical  Report  99-02,  Naval  Medical  Research 
Center,  May  1999. 

5.  P.  K.  Weathersby,  S.  S.  Survanshi,  R.  Y.  Nishi,  and  E.  D.  Thalmann,  Statistically 
Based  Decompression  Tables  VII:  Selection  and  Treatment  of  Primary  Air  and 
N202  Data,  Joint  Technical  Report,  NSMRL  #1182  and  NMRI  92-85,  Naval 
Medical  Research  Institute,  Sept  1992. 

6.  D.  R.  Leitch  and  E.  E.  P.  Barnard,  “Observations  on  No-Stop  and  Repetitive  Air 
and  Oxynitrogen  Diving,”  Undersea  Biomed  Res,  Vol.  9  (June  1982),  pp.  IIS- 
129. 

7.  J.  A.  Hawkins,  C.  W.  Shilling,  and  R.  A.  Hansen,  “A  Suggested  Change  in 
Calculating  Decompression  Tables  for  Diving,”  US  Navy  Med  Bull,  Vol.  33  (1 935), 
pp.  327-338. 

8.  O.  E.  Van  der  Aue,  R.  J.  Kellar,  E.  S.  Brinton,  G.  Barron,  H.  D.  Gilliam,  and  R.  J. 
Jones,  Calculation  and  Testing  of  Decompression  Tables  for  Air  Dives  Employing 
the  Procedure  of  Surface  Decompression  and  the  Use  of  Oxygen,  NEDU  Report 
No.  1  ,Navy  Experimental  Diving  Unit,  1 951 . 

9.  G.  Albano,  “Experimental  Research  on  the  Decompression  of  Man:  I.  Critical 
Safety  Values  of  the  Pressure  Gradient  during  a  Direct  Ascent  to  the  Surface 
without  Stops,”  Med  dello  Sport,  Vol.  1  (1961),  pp.  259-275. 

10.  U.K.  Military  Diving  Manual,  Volume  II,  Change  I,  Book  of  Reference  2806 
(Fareham,  Hants,  U.K:  Commander  in  Chief,  Fleet,  2001). 

1 1 .  G.  R.  Lauckner  and  R.  Y.  Nishi,  Canadian  Forces  Air  Decompression  Tables, 
Report  No.  85-R-03  (Downsview,  Ont.,  Canada:  Defence  and  Civil  Institute  for 
Environmental  Medicine,  1985). 


34 


12.  J.  L.  Meliet,  Les  Tables  de  Plongee  a  I’Air  de  la  Marine  National:  Historique, 

Etude  Statistique,  Propositions  d’Amelioration,  Proces-verbal  No.  03/90  (Toulon 
Naval,  France:  Commission  d’Etudes  Pratiques  d1  Intervention  Sous  la  Mer, 

1990). 

13. H.  D.  Van  Liew  and  E.  T.  Flynn,  A  Simple  Probabilistic  Model  for  Estimating  the 
Risk  of  Standard  Air  Dives.  NEDU  TR  04-42,  Navy  Experimental  Diving  Unit, 

Dec  2004. 

14.  H.  D.  Van  Liew  and  E.  T.  Flynn,  “Direct  Ascent  from  Air  and  N2-02  Saturation 
Dives  in  Humans:  DCS  Risk  and  Evidence  of  a  Threshold,”  Undersea 
Hyperbaric  Medicine,  In  Press. 

15.  S.  S.  Survanshi,  E.  D.  Parker,  E.  D.  Thalmann,  and  P.  K.  Weathersby, 

Statistically  Based  Decompression  Tables  XII:  Repetitive  Decompression  Tables 
for  Air  and  Constant  0.7  ATA  PO2  in  N2  Using  a  Probabilistic  Model,  Naval 
Medical  Research  Institute  Technical  Report  97-36,  Naval  Medical  Research 
Institute,  Aug  1997. 

16.  G.  J.  Duffner  and  H.  H.  Snider,  “Effects  of  Exposing  Men  to  Compressed  Air  and 
Helium-Oxygen  Mixtures  for  12  Hours  at  Pressures  of  2-2.6  Atmospheres,” 

NEDU  Research  Report  1-59,  Navy  Experimental  Diving  Unit,  Jan  1959. 

17.  L.  Wilkinson,  SYSTAT:  The  System  for  Statistics  (Evanston,  IL:  SYSTAT  Inc., 
1990). 

18.  H.  V.  Hempleman,  “History  of  Decompression  Procedures,”  in  P.  B.  Bennett  and 
D.  H.  Elliott,  eds.,  The  Physiology  and  Medicine  of  Diving,  4th  ed.  (Philadelphia, 
PA:  Saunders,  1993),  pp.  342-375. 

19.  J.  Conkin,  V.  Kumar,  M.  P.  Powell,  P.  P.  Foster,  and  J.  M.  Waligora,  “A 
Probabilistic  Model  of  Hypobaric  Decompression  Sickness  Based  on  66 
Chamber  Tests,”  Aviat  Space  Environ  Med,  Vol.  67  (1996),  pp.  176-183. 

20.  H.  H.  Ku,  “Notes  on  the  Use  of  Propagation  of  Error  Formulas,”  J  Res  Natl  Bur 
Stand  Eng  Instr,  Vol.  70C  (October-December  1966),  pp.  263-273. 

21 .  J.  R.  Hays,  B.  L.  Hart,  P.  K.  Weathersby,  S.  S.  Survanshi,  L.  D.  Homer,  and  E.  T. 
Flynn,  Statistically  Based  Decompression  Tables  IV:  Extension  to  Air  and  N202 
Saturation  Diving,  Naval  Medical  Research  Institute  Report  86-51 ,  Naval  Medical 
Research  Institute,  1 986. 

22.  P.  K.  Weathersby,  S.  S.  Survanshi,  and  R.  Y.  Nishi,  “Relative  Decompression 
Risk  of  Dry  and  Wet  Chamber  Air  Dives,”  Undersea  Biomed  Res,  Vol.  17  (July 
1990),  pp.  333-352. 


35 


23.  R.  Ball,  C.  E.  Lehner,  and  E.  C.  Parker,  “Predicting  Risk  of  Decompression 
Sickness  in  Humans  from  Outcomes  in  Sheep,”  JAppI  Physiol,  Vol.  86  (June 
1999),  pp.  1920-1929. 

24.  H.  D.  Van  Liew  and  E.  T.  Flynn,  ’’Use  of  Dive-Outcome  Information  Three  Times 
in  Preparation  of  a  Probabilistic  Decompression  Model,”  Undersea  Hyperbaric 
Med,  Vol.  29  (Summer  2002),  pp.  111-112. 

25.  E.  D.  Thalmann,  Suitability  of  the  USN  MK  15  (VVAL-18)  Decompression 
Algorithm  for  Air  Diving,  NEDU  TR  03-12,  Navy  Experimental  Diving  Unit,  Aug 
2003  (originally  published  as  Final  Report  on  Research  Contract  N0463A-96-M- 
7036,  March  30,  1997). 

26.  E.  D.  Thalmann,  Phase  II  Testing  of  Decompression  Algorithms  for  Use  in  the 
U.S.  Navy  Underwater  Decompression  Computer,  NEDU  Report  1-84,  Navy 
Experimental  Diving  Unit,  Jan  1 984. 

27.  H.  D.  Van  Liew  and  E.  T.  Flynn,  “Probabilistic  Models  for  Standard  Air  Dives: 
Effect  of  Inclusion  of  Saturation  Data,”  Undersea  Hyperbaric  Med,  Vol.  28 
(Supplement  2001 ),  p.  41 . 

28.  C.  D.  Lehner  and  E.  H.  Lanphier,  “Influence  of  Pressure  Profile  on  DCS 
Symptoms,”  in  R.  D.  Vann,  ed.,  The  Physiological  Basis  of  Decompression:  38th 
UHMS  Workshop,  Undersea  and  Hyperbaric  Medical  Society  (Bethesda,  MD: 
UHMS,  1989),  pp.  299-322. 

29.  C.  D.  Lehner,  D.  J.  Hei,  M.  Palta,  E.  N.  Lightfoot,  and  E.  H.  Lanphier, 
“Accelerated  Onset  of  Decompression  Sickness  in  Sheep  after  Short,  Deep 
Dives,”  in  A.  A.  Bove,  A.  J.  Bachrach,  and  L.  J.  Greenbaum,  eds.,  9th 
International  Symposium  on  Underwater  and  Hyperbaric  Physiology  { Bethesda, 
MD:  Undersea  and  Hyperbaric  Medical  Society,  1987),  pp.  197-205. 

30.  R.  D.  Vann  and  E.  D.  Thalmann,  “Decompression  Physiology  and  Practice,”  in  P. 
B.  Bennett  and  D.  H.  Elliott,  eds.,  The  Physiology  and  Medicine  of  Diving,  4th  ed. 
(Philadelphia,  PA:  Saunders,  1993),  pp.  376-432. 

31.  E.  T.  Flynn,  E.  C.  Parker,  and  R.  Ball,  “Risk  of  Decompression  Sickness  in 
Shallow  No-Stop  Air  Diving:  An  Analysis  of  U.S.  Navy  Experience,  1990-1994,” 
in  N.  E.  Smith  and  M.  R.  Collie,  eds.,  Proceedings  of  the  14th  Meeting  of  the 
United  States-Japan  Cooperative  Program  in  Natural  Resources  (UJNR):  Panel 
on  Diving  Physiology,  Panama  City,  Florida,  September  16-17,  1997  (Silver 
Spring,  MD:  U.S.  Dept,  of  Commerce,  NOAA,  National  Undersea  Research 
Program,  [1998]),  pp.  23-38. 


36 


PAGE  INTENTIONALLY  LEFT  BLANK 


CO 

LU 

CO 

< 

O 

CO 

o 

Q 

LL 

O 

CO 

_J 

s 

111 

Q 

>< 

Q 

Z 

UJ 

Q. 

0. 

< 


<0 

Ul 

> 

o 

z 

o 

0) 

CC 

Ul 

0. 

r-~ 
co^r 
O  CD 
cm  w 
— '  co 

L_  -O 

uj  .2 

m  co 

x  Q 

W  c 

Q  .2 

<  CO 
<0 
CD 

i_ 

Q- 

E 

g 
CD 


UJ 

cc 

CL 
U> 

CL 

£o 

> 
> 
cc 


CO 

6 

z  Z 

co 

H  => 

CO  ^ 
UJ  o 
c 
o 

"co 

H— ' 

c 


CO 
< 
o 

CO 

o  © 
O  E 


o 

u. 

o 


3 

o 

o 

T3 

E 

o 


CO  -= 

—I  XJ 

<  i 

uj  05 

q  a 


CD 

to 

c 

o 

c 

CO 

c 

E 

3 


CD 

-o 

O 


CO 

3 

to 

3 

XI 

§ 


£  2 
>>  C 


CO 


w 

CD 


■gl 

5  C  «  II 

‘S  g  >- 
1  -o  co  . . 

cQ-Q 

o  =  E 
£  ©  * 
E  Q-o 
o  ^ 
o  • 


—  CO 
1  © 
CD  o 

i  i-  g 

w  •— 

o 


c 

o 

to 

to 

CD 

i — 

Cl 

E 

§ 

CD 


CD 

Q. 

>* 


_  *  <0 


CO  •— 


c 

CD 


ii  =  E 
'  ©  ©  co 

-+-*  ■*— * 

CO  CO  II 
OO  X 
o  o  ct 


c 

o 

13. 

*L_ 

o 

(/) 

CD 

O 

Bilateral  knee  pain 

Right  knee  pain 

Right  knee  pain 

Left  knee  pain 

Right  knee,  left  ankle  pain 

Right  knee,  left  ankle  pain 

Right  knee  pain,  rash,  mild  nausea,  malaise 

Marked  right  knee  extensor  weakness 

Left  shoulder  pain 

Right  knee  pain 

Mild  knee  pain 

Knee  pain 

Knee  pain 

Ankle  pain 

Left  knee  pain 

Shoulder  pain 

Vertigo,  difficult  breathing,  syncope 

Right  shoulder  pain,  weakness 

Labeled  Type  1;  no  narrative 

Paresthesia,  moderate  fatigue 

Itch,  chest  pain  deep  inspiration 

X 

CC 

> 

> 

>- 

> 

> 

>- 

> 

B 

B 

B 

a 

a 

DCS 

type 

- 

- 

►— i 

i 

>— i 

>— i 

8 

h-H 

- 

- 

a 

>— i 

HH 

i— < 

t-H 

. 

>— h 

Pdcs 

8.3% 

-P 

0s 

CO 

CO 

vP 

O' 

CO 

00 

nP 

o'* 

CO 

CO 

-o 

h- 

| 

vO 

O^ 

q 

CO 

CM 

0s* 

q 

CO 

CM 

nP 

O'' 

CD 

co 

CM 

* 

q 

co 

CM 

^9 

0s 

Is- 

O'' 

CM 

05 

CM 

-9 

q 

CM 

xO 

q 

CM 

-9 

0s 

CD 

T-‘ 

CM 

-9 

o'* 

q 

CM 

-9 

o- 

CD 

CD 

-9 

o'* 

LO 

O 

vP 

O'* 

in 

o 

^0 

O'* 

in 

8.1% 

TDT, 

min 

T— 

1.0 

q 

T“ 

L_ 

1.0 

o 

T— 

o 

o 

o 

T— 

O 

CO 

CO 

05 

cd 

05 

cd 

n 

iq 

CM 

CO 

Bottom 
time,  min 

2,880 

o 

CO 

00 

CM 

o 

CO 

00 

cm" 

o 

txi 

00 

of 

o 

1,440 

1,440 

1,440 

o 

T— 

o 

T — 

o 

■*sJ- 

m 

4,320 

Oi 

N 

“i 

CD 

T“ 

co 

Depth, 

fswq 

25.1 

LO 

CM 

If) 

CM 

1- 

td 

CM 

CD 

CM 

o 

cd 

co 

o 

CO 

CO 

o 

CO 

CO 

O 

CO 

CO 

q 

CO 

co 

CD 

CM 

o 

cd 

co 

to 

OO 

CM 

lo 

O) 

CM 

to 

to 

CM 

to 

05 

CM 

q 

-M- 

CM 

o 

CO 

CM 

o 

co 

CM 

o 

CD 

M- 

94.0 

114.0 

Profile 

number 

CO 

r— 

1 

\ 

1 

CD 

o 

8 

CO 

CM 

i — 

CO 

CM 

Source 

file 

ASATARE 

1 

II 

S 

1, 

1 

ASATDC 

ASATDC 

o 

Q 

H 

< 

CO 

< 

ASATDC 

ASATDC 

ASATNSM 

ASATNSM 

ASATNSM 

ASATNMR 

DC4W 

DC4W 

<: 

O’ 

O 

Q 

D 

o 

D 

z 

to 

CO 

5 

Q 

LU 

EDU1351NL 

Case 

number 

- 

1 

1 

< 


APPENDIX:  DETAILS  OF  DCS  CASES 


CM 

i 

< 


APPENDIX:  DETAILS  OF  DCS  CASES 


,g>  2 


CD 

C  C 
co  to  £ 
Q.  O.  cO 
CD  CD  O- 


I  g>  LffrS 


* 

CD 

c  c 
’co  ‘co 

Q.  CL 
CD  CD  “l. 
CD  CD  CD 
C  C  <D 
^  c 

x:  x: 

tr  cc  _j 


h—  r— 

©  g> 

CD  *- 

c  v_-  £ 

O  CD  S 

x>  -a  S  .£  c 

m  -  3  i_  JO  a 

JB  (0  O  CD  Q-  Q. 

Is-  Q.  SZ  ~0  a  -> 

s  0  w  3  e  | 

©  2  .&  ?  I  i 

d  e/>  £  cd 

C 

co  "aj  "a3  cd 

Q _ I  _J  _1 


gLri-gLrcLgL--*  F  ^  -1  Q-  -1  -J  H 
>zzzz  >->->>>->->>- 


*  *  *  *' 


CM  CM  CM  CM  CM 

C/3  C/3  C/3  C/3  C/3 

03  CD  CD  03  03 

•c}-  -cf  ^3- 
CO  00  00  00  CO 

3  3  3  3  3 

Q  Q  O  Q  3 

UJ  ua  UJ  LU  LU 


N  CO  O)  O  i- 
ID  LO  LO  CD  CO 


-3 

a  •- 

— 

1— H 

►— < 

k— <  >-H 

Sp 

>P 

nO 

vP  -sO 
O'  0s 

LO  LO 

LO 

LT5 

m  in 

o  o 

o 

d 

o  o 

co  co 

co 

co  co 

co  co 

CD 

CO 

CD  CD 

T“ 

T~ 

T-  T- 

o  o 

o 

o 

o  o 

CM  CM 

CM 

CM 

CM  CM 

r-  h- 

r- 

r- 

N  h- 

o  o 

o 

o 

O  O 

o  o 

o 

o 

d  o 

^  ■M' 

O  1- 

C\J 

CO 

■M-  in 

CO  CO 

CO 

CO 

CO  CO 

CM  CM 

CM 

CM 

CM  CM 

C/3  C/3 

C/3 

C/3 

C/3  C/3 

0)  05 

03 

CD 

03  03 

-cT  ■M- 

CO  CO 

CO 

CO 

00  CO 

3  3 

3 

3 

3  3 

D  D 

D 

Q 

Q  Q 

UJ  UJ 

UJ 

UJ 

UJ  UJ 

CM  CO 

LO 

CD 

CO  CO 

CD 

CO 

CO  CO 

O  O 


Ai  if*  f n  rs^  cn  (T5  O  t—  C\i  CO  "M"  LO  CO  Is--  co  05  O  t-  C\J  CO  ^  LO  CD  CO  05 

(£)0(DCDCDCDCOCDNNNNNr^NNNNCDCOCOCO  CO  CO  CO  CO  00  00 


APPENDIX:  DETAILS  OF  DCS  CASES 


A-4 


